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DIRECTIONAL  UNDERWATER  NOISE  ESTIMATES  - 
THE  DUNES  MODEL 


R.W.  Bannister,  D.J.  Kewley  and  A.S.  Burgess 


SUMMARY(U) 

The  modelling  of  the  directionality  of  underwater  noise  due  to  wind -generated 
noise  sources  and  shipping  is  presented.  The  formulations  are  incorporated  in 
the  Directional  Underwater  Noise  Estimates  -  DUNES  model.  It  provides 
estimates  of  omnidirectional,  vertical,  horizontal  and  three  dimensional 
directional  noise  versus  frequency.  The  model  includes  features  hitherto 
unknown  in  models  of  this  type  ie  high  latitude  and  coastal  slope  enhanced 
wind  noise.  The  model  emphasises  the  calculation  of  noise  due  to  the  natural 
environment  and  therefore  shipping  contributions  have  to  be  entered 
explicitly,  unlike  some  other  models  which  incorporate  extensive  shipping 
databases.  Long  range  wind  noise  is  considered  to  be  described  by  storms  over 
a  finite  area  and  bearing  or  along  continental  shelves.  Shipping  is  described  by 
individual  ships  or  shipping  lanes  over  continental  shelves  or  deep  ocean.  The 
model  has  been  used  on  a  range  of  computers  such  as  the  IBM  PC,  VAX  1 1/750 
and  IBM  3033.  Examples  of  use  for  Version  2.3  are  p  ‘  ""led. 
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1.  INTRODUCTION 

The  waters  surrounding  Australia  and  New  Zealand  are  relatively  free  from  ship 
traffic.  As  a  result,  natural  noise  sources  assume  significance  in  the 
frequency  band  20  to  200  Hz.  This  is  in  contrast  to  many  areas  in  the 
northern  hemisphere  where  shipping  noise  dominates  this  part  of  the  spectrum. 
In  spite  of  these  differences,  however,  there  are  many  features  of  the 
measured  three-dimensional  ambient  noise  field  which  appear  similar  in  both 
hemispheres.  For  example  there  is  a  broad  spectral  peak  in  the  vicinity  of 
50  Hz(ref.l)  and  the  maximum  of  omnidirectional  noise  occurs  around  SOFAR 
channel  depths (ref . 2 , 3) .  In  the  vertical  plane,  there  is  a  peak  in  energy 
received  at  angles  near  to  the  horizontal (ref . 2 , 3)  and  horizontal  noise 
directionality  shows  considerable  azimuthal  dependence. 

In  the  northern  hemisphere,  these  attributes  have  been  shown  to  be  dominated 
by  the  presence  of  shipping,  often  enhanced  by  low- loss  transmission  paths 
which  occur  when  vessels  are  at  high  latitudes  or  over  sloping 
bathymetry  (. ref .  4 )  .  Similar  characteristics  which  are  measured  in  southern 
oceans  (.  ref .  3 , 5 , 6 )  do  not  easily  fit  the  same  explanation.  There  are  fewer 
ships  overall,  negligible  numbers  suitably  placed  above  bathymetric  slopes  and 
practically  none  at  all  at  high  latitude. 

The  apparent  inconsistency  between  interpretation  of  the  broad  features  of 
northern  and  southern  hemisphere  noise  data  has  challenged  us  to  attempt  to 
further  understand  the  role  and  characteristics  of  wind  generated  ambient 
noise.  Since  the  difficulties  appear  to  be  fundamental  it  was  decided  to 
construct  an  ambient  noise  model  which  emphasises  in  simple  form  the  various 
physical  processes  in  operation.  The  wind  generated  noise  source  is  regarded 
as  being  distributed  over  the  ocean  surface  with  simple  latitude  and  seasonal 
variation.  The  model  simulates  the  acoustic  'view'  at  a  receiving  point  and 
traces  various  transmission  processes  back  to  an  interaction  with  the  sea 
surface  source.  Each  transmission  family  his  characteristic  physical 
properties  and  can  be  associated  with  a  specific  set  of  vertical  and 
horizontal  angles.  This  model,  known  as  DUNES  (Directional  Underwater  Noise 
Estimates  or  more  colloquially  Down  Under  Noise  Estimates)  has  been 
benchmarked  against  a  set  of  vertical  noise  measurements  from  southern  ocean 
environments . 

The  model  has  been,  and  is  continuing  to  be  modified  in  response  to  revised 
interpretations  of  noise  mechanisms  and  paths.  Consequently  this  report 
attempts  to  describe  Version  2.3  which,  as  seen  in  the  text,  can  be  further 
improved  by  removing  assumptions  and  adding  more  features.  Therefore  by  the 
time  this  Technical  Note  appears,  new  research  may  require  further  changes. 

This  report  is  organised  into  two  main  sections.  Part  one  describes  (and 
derives  where  necessary)  the  algorithms  used  in  the  noise  model.  Part  two  is 
a  users  guide  for  the  software  implementation  of  the  model  and  describes  how 
to  establish  and  arrange  parameters.  Additional  notes  and  Worksheet  are  given 
in  the  Appendices. 


2.  THE  DUNES  MODEL  ALGORITHMS 


2 . 1  Overview 

The  DUNES  model  estimates  the  noise  arriving  at  a  receiver  as  the  sum  of 
contributions  associated  with  different  transmission  types  as  indicated  in 
figure  1.  Each  transmission  type  can  be  assigned  a  vertical  angle  interval 
at  the  receiver  as  indicated  in  figure  2  (derived  from  reference  7).  The 
angle  defining  the  change  from  the  direct  to  BB  components  is  represented 
in  figure  2  by  the  critical  angle.  Rays  arriving  at  angles  greater  than 
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this  critical  angle  undergo  much  larger  losses  upon  reflection  at  the  sea 
bottom  than  for  smaller  angles.  Its  value  can  be  found  from  the  angular 
dependence  of  bottom  loss,  see  Appendix  III.  The  DL'N'ES  model  also 
allocates  horizontal  sections  for  the  various  components.  It  can  be 
assumed  that  local  (or  direct),  RSR  and  BB  components  arrive  uniformly  from 
all  horizontal  directions  (unless  major  topographic  blockage  occurs  within 
100  km  of  the  receiver).  The  RR  component  how’ever,  arises  from  distant 
noise  sources  either  over  seamounts,  over  continental  dopes  or  from  high 
latitude  and  the  azimuth  can  be  determined  from  geometry. 

The  following  sections  will  develop  and  discuss  the  algorithms  used  in 
DUNES  to  predict  the  noise  components  illustrated  in  figures  1  and  2.  The 
values  of  the  vertical  angle  assignments  shown  in  figure  2  correspond  to 
the  case  of  deep  water  with  a  depth  excess  in  the  sound  speed  profile.  For 
cases  of  bottom  limited  or  no  deep  sound  channel  propagation  these 
assignments  will  change  and  whole  intervals  can  be  eliminated.  The  values 
are  discussed  in  Sections  2  and  3. 

2.2  Source  ievel  of  wind  noise 

2.2.1  Source  spectrum 

Wind  generated  noise  arises  from  a  distributed  source  near  to  the  ocean 
surface.  The  source  level  is  described  quantitatively  in  decibels  re 
1  yPa2/Hz/m2  of  surface  area  at  one  metre.  This  has  been  estimated 
recently  by  Burgess  and  Kewley ( ref  .  5 )  and  by  Wilson ( ref . 8 )  .  Wilson's 
source  levels  are  derived  from  omnidirectional  noise  measurements 
selected  from  locations  where  direct  path  conditions  dominate.  Burgess 
and  Kewley,  have  derived  an  alternative  technique  to  directly  measure 
the  surface  source  level  using  upward  and  downward  steered  beams  from  a 
vertical  array.  This  method  specifically  accounts  for  the  effects  of 
bottom  ioss  on  the  measured  vertical  noise  levels.  The  Burgess  and 
Kewley  approach  estimates  effective  downwards  monopole  levels  ie  the 
product  of  the  source  level  and  its  downwards  directivity  function. 
This  is  important  in  the  present  study,  since  surface  dipole  effects  are 
explicitly  included  in  the  propagation  loss  models  and  so  effective 
monopole  source  levels  must  be  used. 

The  analytic  form  of  source  level  used  in  the  model  is  obtained  from  an 
updated  version  of  Burgess  and  Kewley 's  curves (ref . 9 )  which  incorporated 
additional  data  obtained  from  the  Indian  Ocean(ref . 10) .  The  expression 
for  the  source  level  at  steep  angles  used  here  is 


S(dB  re  1  yPa2/Hz  -  m2  at  1  m) 


=  St  -  2.8(S2  -  Sl)  +  l.o6(S2  -  Sjjlog  f 


U) 


where  Sj 

S2 

f 

w 


37  +  12  log  W 

25  +  20  log  W 

frequency  (Hz) 
wind  speed  (kn) 


In  addition,  values  of  S  below  50  Hz  and  between  200  Hz  and  1  kHz  are 
held  constant.  The  resulting  curves  are  shown  in  figure  3(a) 
superimposed  on  the  source  level  measurements  of  Burgess  and 
Kewley ( ref . 9 ) .  (It  is  shown  in  Appendix  V  that  the  error  in  converting 
from  units  of  dB  re  1  yPa2/(Hz  -  sr)  to  dB  re  1  yPa2/(Hz  -  m2 )  at  1  m  is 


3 


VSRL-TN-34/89 


less  than  0.3  dB  so  both  sets  of  units  are  virtually  interchangeable). 
At  higher  frequencies  the  values  above  1  kHz  are  reduced  by 
11  dB/decade.  This  reduction  leads  to  the  observed  17  dB/decade  falloff 
of  omnidirectional  noise  level.  At  very  low  frequencies  the  data 
summary’  of  Kibblewhite  and  Ewans  ( ref .  1 1 )  was  used  to  estimate  an 
additional  noise  term 


S  =  16.6  -  50  log(f/50)  +  10  log(W/10) 


to  be  combined  with  S  in  equation  (1)  by  power  summation  to  give  values 
of  source  level  down  to  10  Hz.  It  is  interesting  to  note  a  major 
difference  between  equation  (1)  and  the  results  quoted  by  Wilson! ref . o j . 
Wilson's  source  level  curves  have  a  uniform  wind  speed  dependence 
(45  log  (wind  speed))  over  the  whole  frequency  band  (10  to  1000  Hz). 
However,  the  expression  derived  from  Burgess  and  Kewlev,  has  a  12  log 
Iwind  speed)  variation  below  50  Hz  and  20  log  (wind  speed;  variation 
above  200  Hz.  The  differences  between  these  two  estimates  can  be 
resolved  if  the  effects  of  source  depth  are  included.  This  point  is 
discussed  in  the  next  section. 

2.2.2  Wind  noise  source  depth 

The  source  level  estimate  made  by  Burgess  and  Kewlev  is  based  on  a 
measurement  using  nearly  vertical  beams.  Thus,  horizontally  arriving 
energy  is  specifically  eliminated  and  results  are  free  from  major 
surface  dipole  loss  effects.  The  Burgess  and  Kewlev  source  level  is 
therefore  an  effective  monopole  value  and  so  surface  dipole  effects  at 
low  angles  must  be  included  in  the  long  range  transmission  loss  estimate 
used  to  assess  received  noise.  Indeed,  it  will  be  shown  iater  that  for 
this  reason  the  received  level  is  dependent  upon  the  monopole  source 
depth.  Consequently,  it  is  necessary  to  estimate  the  source  depth  of 
wind  generated  noise. 

The  most  desirable  approach  to  this  problem  would  be  to  start  with  an 
understanding  of  the  physical  mechanism  of  wind  generated  noise  and  from 
this  to  estimate  an  effective  source  depth.  Unfortunately  the 
mechanisms  are  not  clearly  established.  This  has  left  us  with  the  need 
to  find  an  empirical  solution  by  comparing  DUNES  predictions  with 
measured  data,  particularly  at  frequencies  less  than  100  Hz.  This 
procedure  has  indicated  clearly  that  the  source  depth  must  be  both  wind 
speed  and  frequency  dependent.  Effective  source  depths  inferred  in  this 
manner  are  shown  in  figure  3(b). 

Some  rationale  can  be  applied  to  these  results  in  physical  terms.  The 
source  model  implies  a  process  which  extends  deeper  with  increasing  wind 
speed.  It  is  reasonable  to  expect  both  turbulence  and  entrained  bubble 
activity  to  extend  to  deeper  depths  under  conditions  of  high  surface.' 
wave  activity  than  under  calm  conditions.  A  possible  mechanism  has  been 
suggested  by  Carey  and  Browning) ref .  12) .  Turbulent  pressure 
fluctuations  which  are  coupled  to  the  water  from  air  turbulence,  may 
oscillate  entrained  bubbles.  The  bubbles  then  become  sound  radiators. 
Since  bubble  clouds  are  observed  in  the  ocean(ref . 13)  at  depths  which 
depend  on  wind  speed,  this  may  well  be  a  relevant  mechanism.  The  source 
consists  of  a  series  of  incoherent  line  radiators  oriented  vertically 
whose  length  changes  with  time  but  varies  consistently  with  wind  speed. 
The  time  variability  is  reduced  by  geographic  averaging  over  a  large 
surface  area  of  the  ocean. 
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Using  the  empirical  relationship  between  significant  wave  height  and 
wind  speed  as(ref.l4) 


V,  =  0.00133  W  2  ' 5 

h  s 

where  =  significant  wave  height  in  metres 

W  =  wind  speed  in  knots, 

a  source  depth  (hg)  of  five  times  the  significant  wave  height  is  found 
to  produce  DUN:ES  predictions  consistent  with  measurement.  Thus 


h  =  0.00665  W  2 ' 5 

s  s 


( 1  ( a ) ) 


This  expression,  perhaps  coincidentally,  also  fits  the  spread  oi  data  n 
reference  13  of  bubble  depth  versus  wind  speed  up  to  20  kn . 

However,  there  is  an  additional  factor  which  is  related  to  the  surface 
interference  effect.  It  can  be  shown  that  the  first  interference 
maximum  as  a  function  of  depth  occurs  at  a  source  deptn  of 
hg  =  0.25/sin  0^  (wave  lengths)  (where  0q  is  the  surface  grazing  angle). 

It  is  assumed  that  once  the  source  depth  has  reached  this  limit,  further 
depth  increase  will  not  significantly  affect  the  radiated  energy.  Thu 
the  algorithm  used  in  the  DUNES  model  is  for  the  source  depth  to  vary  as 
in  equation  (1(a))  with  a  lower  limit  of  X/4.  Figure  3(b)  shows  this 
graphically. 

The  assumption  of  X/4  as  the  limit  for  source  depth  fulfills  the  authors 
expectation  that  the  source  of  noise  w ’ 1 1  be  at  shallower  depths  as  the 
frequency  increases.  However  at  the  higher  frequencies  considered  in 
the  model  (up  to  10  kHz)  this  depth  is  comparable  to  the  variation  of 
the  sea  surface  away  from  the  horizontal.  The  assumption  of  a  pure 
dipole  surface  interference  effect  with  a  flat  surface  will  therefore  be 
unlikely  to  be  valid  for  most  wind  speeds.  However,  as  will  be  seen 
later,  the  dipole  does  not  contribute  to  the  longer  range  noise  at  these 
frequencies  and  an  effective  monopole  source  dominates  most  of  the  noise 
levels  at  low  receive  angles  IBB  and  RSR).  Thus  the  X,/4  is  shown  to  be 
not  a  critical  factor  at  the  higher  frequencies.  As  stated  before, 
physical  understanding  of  the  sources  of  noise  will  eventually  tell  us 
whether  these  assumptions  are  valid.  Thus  a  plane  of  sources  at  X/4  may 
be  too  simple  an  approximat ion  to  represent  a  volume  of  sources  over  a 
larger  depth  range. 

2.3  The  refracted  -  refracted  (RR1  components 

The  existence  and  importance  of  the  RR  components  has  been  known  for  some 
time  and  stressed  in  recent  papers ( re f . 2 ,4 , 15 ) .  The  RR  component  requires 
horizontal  gradients  in  either  bathymetry  or  sound  speed  to  carry  energy 
from  surface  sources  to  the  deep  RR  duct.  These  conditions  are  met  when 
noise  sources  (either  shipping(ref . 2)  or  wind(ref,16)  are  over  bathymetric 
slopes  or  at  high  latitudes  where  the  SOFAR  channel  axis  approaches  the 
surface (ref,  15  ) )  . 
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Situations  of  this  sort  are  conveniently  approached  by  integrating  the 

noise  contribution  of  a  rectangular  distributed  source  element  over 

appropriate  range  intervals.  As  discussed  in  reference  16  the  cylindrical 

transmission  loss  ( L  )  of  the  form 

KK 

Ldd  =  70+10  log  R  +  aR  (2) 

KK 

where  R  =  Range  (km) 

a  =  Water  column  attenuation  appropriate  to  the  sound 

source  (dB/km) 

can  be  assumed. 

References  16  and  17  show  that  the  received  noise  level  VI  )  due  to  this 

KK 

component  is  estimated  by 

1DD  =  S„  +  10  log  W  -  70  +  3  +  K  -  aR/xi/4R2  +  1  13) 

KK  K 

10  log[ln(xi/R  +  /xj/R2  +  1)] 

half  the  length  of  the  symmetrical  noise  lane  fkm) 

the  perpendicular  range  of  the  noise  lane  from  the 

receiver  (km) 

the  width  the  noise  lane  fkm) 
the  source  level  per  kilometres2 

S  +  60 

Note  that  the  noise  lanes  used  in  DUNES  are  assumed  to  be  symmetrical  about 
the  point  of  closest  approach  and  therefore  equation  (3)  differs  from 
equation  (A4)  of  reference  16  by  3  dB . 

Equation  (3)  can  be  applied  to  distributed  wind  noise  sources  at  high 

latitude  or  over  bathymetric  slopes  as  discussed  in  reference  16.  The  same 
equation  applies  to  shipping  in  similar  situations  if  the  source  level  of 
shipping  is  assumed  to  be  uniformly  distributed  over  a  ship  lane  -  this  is 
discussed  in  Section  2.9.2  The  RR  contribution  is  received  over  vertical 

angles  of  ±0  =  cos  1  (C  /C  )  where  C  ,  C  are  sound  speeds  at  the  receiver 

o  o  s  o  s 

and  surface.  DUNES  does  not  allow  an  RR  contribution  unless  C  >  C  and 

s  o 


High - lat itude ,  wind-generated  noise  has  a  seasonal  and  latitude  dependence 
as  discussed  in  reference  16.  The  frequency  spectrum  of  this  component 
f  figure  11)  is  such  as  to  confine  its  significance  to  the  total  noise  to 
frequencies  below  200  Hz.  This  effect  is  due  to  the  water  column 
attenuation  of  the  high  frequencies  over  the  very  long  ranges  to  the  high 
latitude  sources.  Note  that  the  high  latitude  wind  calculations  in  DUNES 
use  the  source  levels  presented  in  reference  16.  These  are  similar  to 
equation  (1)  at  10  kn  but  have  greater  values  at  higher  wind  speeds.  These 
values  were  used  so  that  DUNES  would  be  consistent  with  the  earlier 

publication.  To  take  into  account  surface  decoupling,  an  f  1  loss  is 
applied  to  I  ^  below  40  Hz.  As  discussed  in  later  sections  the  loss  term 
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depends  upon  source  directionality  and  depth,  giving  rise  to  an  increasing 
loss  for  decreasing  frequency.  The  loss  term  used  is  a  first 
approximation . 

The  additional  contribution  from  sources  over  bathymetric  slopes  extends  to 
higher  frequencies  C figure  11).  This  is  a  consequence  of  the  shorter 
average  distance  between  source  and  receiver  producing  less  high  frequency 
attenuation  in  the  water  column.  Thus,  the  slope  enhanced  component  is 
likely  to  dominate  in  the  RR  angles  above  200  Hz.  This  component  is,  of 
course,  very  dependent  on  the  location  of  the  receiver  relative  to 
bathymetric  features. 

2.4  The  RSR  component  -  transmission  loss 
2.4.1  General  law 

The  RSR  component  does  not  have  any  bottom  interactions  so  it  can  be 

reasonably  expected  to  have  a  cylindrical  spreading  loss  (LDCD)  plus 

KoK 

water  column  attenuation, 
viz 


Ldc_  =  70+10  log  R  +  aR  (4) 

KoK 

where  the  value  of  a  is  appropriate  to  the  receiver  location.  The 
constant  70  is  the  generally  accepted  empirical  constant  obtained  for 
RSR  paths  to  deep  sources  and  receivers.  In  the  present  case,  however, 
the  effective  source  depth  is  shallow  and  surface  dipole  effects  need  to 
be  included. 

This  is  done  by  integrating  the  surface  dipole  response  over  the  angle 

limits  relevant  to  RSR.  These  are.  between  zero  and  0,  =  cos  1  (C  /C„). 

o  B 

The  latter  angle  describes  the  launch  angle  of  a  ray  at  the  receiver 
depth  which  grazes  the  ocean  bottom.  C  is  the  water  speed  of  sopnd  at 

D 

the  bottom.  The  dipole  intensity  (I  )  as  a  function  of  surface  angle  0 
and  source  depth  z  (wavelengths)  is 

Iz  =  4| P| 2  sin2  X 

where  | P | 2  =  monopole  power  level 

and  X  =  phase  difference  between  the  direct  path  ray  and  surface 

reflected  path  ray 

Assuming  isospeed- conditions  at  the  source  depth: 


X  =  2ir  z  sin  0 

-  2ttz0  (if  sin  0=0) 


Thus  the  integrated  surface  dipole  loss  is 


7 


W'SRL-TN-34/89 


Dipole 


sin2  (2tt  z  0)  dQ 


=  2f);  -  sin  4n  z  01/(2n  z) 


When  this  is  compared  to  the  equivalent  monopole  energy  loss  over  the 
same  angle  limits  (namely  0!)  for  which  the  simple  loss  law 

(equation  (4))  applies,  then  the  resulting  decoupling  factor  relative  to 
equation  (4)  conditions  is  given  by 


drsr  1dipole/1monopole 


=  (20!  -  sin  4it  z  0!/(2tt  z))/Qx 


(5) 


=  2  -  sin(4ir  z  0x)/(2ir  z  0 x ) 


This  is  graphed  in  figures  4  and  5,  from  which  the  following  points  can 
be  made. 

2.4.2  The  effect  of  source  depth 

The  process  is  very  dependent  upon  source  depth.  If  the  source  depth  of 
wind  generated  noise  is  assumed  constant  (in  metres)  then  there  is  a 
6  dB/octave  low-frequency  roll-off.  For  example,  assuming  a  7.5  m  depth 
this  gives  typical  losses  of  -2  dB  at  200  Hz;  -8  dB  at  100  Hz,  -14  dB  at 
50  Hz  and  -20  dB  at  25  Hz  (figure  5). 

Unfortunately,  the  experimental  sites  used  to  benchmark  the  model  do  not 
exhibit  clear  evidence  of  the  RSR  component.  It  has  thus  not  been 
possible  to  check  the  above  formulation  directly.  At  present  the 
assumed  source  depth  is  found  by  making  it  wind  speed  and  frequency 
dependent  as  discussed  in  Section  2.2.2  and  figure  3(b). 

2.4.3  The  effect  of  the  dipole  null 

As  the  surface  angle  approaches  zero,  the  response  of  the  surface  dipole 
becomes  very  low.  When  the  local  sound  speed  profile  is  such  that  RSR 
rays  are  confined  to  near  zero  angles,  the  contribution  is  significantly- 
reduced.  This  is  clear  in  figures  4  and  5.  The  dipole  null,  however, 
is  not  likely  to  be  perfect  and  some  lower  limit  to  the  null  should  be 
postulated  on  physical  grounds.  Null-filling  will  be  caused  by 
scattering  and  the  rays  responding  to  the  detailed  shape  of  the 
near-surface  sound  speed  profile.  In  the  absence  of  any  properly 
developed  theory  an  arbitrary  lower  limit  of  -10  dB  is  set  by  adding  0.1 
to  the  decoupling  factor.  This  decision  is  not  robust  since  the 
magnitude  and  frequency  content  of  the  RSR  component  is  sensitive  to 
this  procedure.  Relevant  data  are  not  available  to  check  results. 
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2.4.4  The  attenuation  coefficient 

The  value  of  a  in  equations  12)  and  (4)  represents  all  relevant  range 

dependent  losses  in  excess  of  cylindrical  spreading.  Clearly,  water 

column  attenuation  (a  )  is  involved,  which  in  the  South  Pacific! ref . 18) 

w 

is  given  by 


a  =  0.07  f  */(l  +  f  *)  +  0.011  f  *  dB/km  (6) 

W  K  K  K 

where  f^.  =  frequency  in  kilohertz.  The  water  column  attenuation  model 

used  in  DUNES  depends  upon  ocean  with  the  expressions  given  in 
reference  19. 

A  further  frequency  independent  term  is  frequently  added  to  account  for 
scattering  which  directs  energy  out  of  water  column  modes  into  ,  the 
bottom.  According  to  Guthrie's  mode  1 ( ref . 20 j  the  oceanographic 
scattering  component  has  the  form 


a  -  1 .  1  x  10  1  e  X/^  dB/km 

os 


(“ ) 


where  h  is  the  characteristic  depth  (200  m)  and  X  is  the  depth  Cm)  on  a 
propagation  path.  This  is  an  empirical  result  based  on  analysis  of  bomb 
arrival  time  data.  Thus  it  can  be  seen  that  the  attenuation  of  ray 
paths  which  traverse  near  to  the  surface  will  be  higher  than  those  which 
are  confined  to  deeper  parts  of  the  water  column.  This  process  can  be 
expected  to  scatter  energy  into  higher  angles  -  some  of  which  being  lost 
into  the  bottom  (hence  the  extra  attenuation  in  water  borne  components  i  . 
Not  all  this  energy  is,  of  course,  lost  to  the  sea  floor.  However,  it 
does  imply  a  transfer  of  energy  into  bottom-bounce  modes. 

The  investigation  of  these  effects  will  require  a  rather  specific 
experiment  and  is  beyond  the  precision  available  with  data  available  to 
the  authors.  Although  these  possibilities  are  recognised,  the  model 
does  not  at  present  attempt  to  emulate  them.  A  value  of  0.0009  dB/km  is 
added  to  the  water  column  attenuation  for  all  oceans. 

2.5  The  RSR  component  -  noise  contribution 

Gathering  together  equations  (4)  and  (5)  above,  the  transmission  loss  law 

for  RSR  including  dipole  and  monopole  contributions  is  ( in  decibels  terms) 


Ldcd  =  70+10  log  R  +  aR  -  10  log  D 

KoK 


where  D 


Dipole  decoupling  term  modified 
2.1  -  sin(4ir  z  0x)/(2tt  z  0x) 


In  intensity  terms,  the  loss  function  is: 


!rSr  De  /Car) 


(8) 


(9) 
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where  10  log  a  =  70,  r  =  R 

and  a'  =  attenuation  in  nepers /ki lometres . 

Consider  an  element  of  the  ocean  surface  as  in  figure  6.  The  area  of  the 
annular  element  (dAE)  is 

dAE  =  2Trr  dr  (10) 

Hence  the  noise  contribution  of  this  element  at  the  receiver  is  (using 
equations  (9)  and  (10)) 

diRSR  =  D2ir/(a)  s^  exp(-a'r)  dr 

where  sR  =  surface  source  level  per  kilometres2. 

Hence  the  total  received  noise 

r  i 

i'RSR  =  D2tt/  ( a )  sR  ^  exp(-a'r  )  dr 

r 

o 

(11) 

=  D2ir/(aa')  [expt'-a'r^  -  expi-a'rj)] 

where  rQ  and  rj  are  the  initial  and  final  ranges  of  integration.  In  decibel 
units 

IrSR  =  \  +  8  *  A  -  10  log  a'  +  10  log  D 

+  10  log(exp(-a '  r  ) -exp( '-a*  r  x ) ,)  (12) 

where  SR  =  10  log  sR,  8  =  10  log  (2tt)  and  A  =  10  log  a.  (Note  a1  =  0.23  a, 

where  a  is  in  decibels/kilometres.)  The  characteristics  of  the  expression 
can  be  judged  by  reference  to  figure  7.  Here  IDCD  is  plotted  as  a  function 

of  rj  (ignoring  the  dipole  term  10  log  D)  and  shows  the  sensitivity  of  the 

RSR  prediction  to  integration  range.  At  300  Hz  the  prediction  has  got  to 
within  3  dB  of  its  asymptotic  value  (rj  =  °° )  at  a  range  rj  =  200  km.  At 

50  Hz  the  3  dB  point  does  not  occur  until  a  range  =  4000  km.  Hence,  the 

low  frequency  noise  contribution  due  to  RSR  extends  to  considerable  ranges 
from  the  receiver.  The  model  implementation  allows  for  bathymetric 
blockage  to  limit  this  range  in  each  of  four  quadrants  to  account  for 
practical  geographical  effects. 
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In  terms  of  sensitivity  to  the  initial  range  of  integration  (rQ),  figure  7 

shows  this  to  be  small  for  r  less  than  200  km.  The  model  therefore  uses 

o 

the  arbitrary  initial  range  of  20  km  to  represent  the  beginning  of  RSR 
noise  contribution.  At  frequencies  above  1  kHz  the  sensitivity  will  be 
greater  but  the  total  contribution  will  be  much  lower  than  that  at  lower 
frequencies  and  therefore  the  errors  with  respect  to  the  total  noise  are 
not  important. 

The  addition  cf  the  dipole  term  in  equation  (12)  reduces  levels  shown  in 
figure  7  by  typically  10  dB .  This  dipole  term  is  most  sensitive  to  the 
source  depth  parameter  (z).  As  discussed  earlier,  it  has  been  found  that 
source  depth  needs  to  be  both  wind  speed  and  frequency  dependent  before 
predicted  values  at  low  frequencies  agree  with  measurement.  At  high 
frequencies  the  monopole  contribution  is  found  to  dominate  and  z  is  of  less 
importance . 

The  resultant  RSR  noise  characteristics  are  summarised  in  figure  11  for 
various  wind  speeds  and  frequencies.  The  energy  is  received  over  the  small 

vertical  angle  band  cos  1  (C  /C  i  to  cos  1  (C  /CL)  where  C  ,  C  and  CL  are 
°  o'  s'  oB  os  B 

sound  speeds  at  the  receiver,  surface  and  bottom  respectively  (as  in 

figure  2).  The  contribution  is  isotropic  unless  bathymetric  effects  limit 

the  propagation  from  certain  azimuths.  There  is  no  RSR  noise  unless 

Cs<CB- 

2.6  The  bottom-bounce  component  -  transmission  loss 

In  the  transmission  region  where  acoustic  mode  stripping  occurs  due  to 
bottom  interaction  (Weston( ref . 21 ) ,  Denham(ref . 22) ) ,  propagation  loss  is  of 
the  form 


L_d  =  A  +  15  iog  R  +  aR 
bb 


(13) 


However,  where  the  source  is  near  to  the  surface,  transmission  loss 
approaches  a  25  log  R  law  due  to  the  influence  of  the  surface  dipole. 
Denham(ref . 22)  has  estimated  that  transmission  loss  for  the  bottom  bounce 
component  for  near  surface  sound  sources  in  the  mode  stripping  region  is 
(approximately) 


lgB  =  1/a  r"3'2  .  gQ(r)  .  exp(-a’r) 


(14) 


where  g  fr') 
k 

o 

h 

s 

H 

Q 

a 

s 


i 

a 


1  -  exp{-(kohs)2  H/(rQ)} 

2tt/X  (m  l) 
source  depth  (m) 

water  depth  (km) 

a  /20  log  e 
s 

57  x  0.05  (f/25) 1  ’ 5  =  slope  of  bot .  loss  (dB/rad^) 
3  (f/25)1'5 

water  attenuation  (nepers/km) 


Denham  showed  reliable  predictions  for  experimental  data.  This 
transmission  loss  law  is  also  found  to  predict  the  results  of  PE 
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calculations  reasonably  well  over  frequencies  and  distances  of  interest 
(figure  8).  The  calculations  were  made  for  a  4748  m  deep  bottom  limited 
propagation  case  with  a  realistic  geoacoustic  model  for  the  sea  floor.  The 
sound  speed  profile  was  from  the  Indian  Ocean  with  C  -  CD  =  9.5  m/s.  Thus 

a  realistic  sound  speed  profile  case  also  can  be  approximated  by 
equation  (14).  The  bottom  loss  slope  parameter  (a^)  has  been  found  to  be 

of  prime  importance.  The  simple  form  used  (as  above)  attempts  to  distil 
global  frequency  and  geographic  variability  into  one  simple  expression.  An 
expression  for  ag  can  be  der ived(ref . 23 , 24 )  giving 


as  =  cfn  (dB/rad) 


where  c  is  a  variable  using  density  and  sound  speed  ratio  terms  and  the  sea 
floor  sediment  compressional  wave  attenuation  coefficient.  The  value  of  n 
for  shallow  water  sediments  is  usually  between  0  and  1 1 ref . 23 , 24) .  We  have 
chosen  n  =  1.5  based  upon  examination  of  bottom  loss  data  for  deep  oceans. 
Final  usage  of  the  noise  model  has  shown  that  regional  changes  in  this 
parameter  are  justified  although  no  scheme  has  currently  been  devisea  to 
quantify  this.  The  current  model  allows  the  user  to  select  a  constant 
multiple  of  the  above  simple  expression  for  ag  ana  retains  tne  same 

frequency  law  at  all  locations. 

To  see  the  essential  loss  law  operating,  the  following  simplification  is 
made (ref . 22) : 


g  (r)  =  (k  h  ') 2  H/(rQl 

o  os  • 

when  r  is  greater  than  6.8  H/OgO2  and  0^  is  the  critical  angle  on  the 
bottom  loss  profile. 

Hence 

IgB  =  H(k0  hs)2/(Qa)  r  5/2  exp(-a’r)  (15) 

or  in  decibels 

Lgg  =  A  +  25  log  r  +  ar  -  10  log  (4n2  h2  H/Q) 

where  h  is  the  source  depth  in  units  of  wavelength. 

Assuming  A  =  63 


Ldd  =  37+25  log  r  -  10  log  (h2  H/a  )  +  ar 

DD  S 


(Note  a'  =  0.23a  where  a  is  in  decibels/kilometres) 

Thus  we  can  see  that  for  near  surface  sound  sources  the  propagation  loss 
has  a  25  log  R  law  with  a  dependence  on  source  and  water  depth. 
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2.7  The  bottom  bounce  component  -  noise  contribution 
2.7.1  Dipole  contribution 

The  transmission  loss  for  the  bottom  bounce  component  including  surface 

dipole  effects  (equation  (16))  is  in  the  form  • 


Loss  =  kr  5/2  exp(-a'r) 


(17) 


where 


k  =  4tt2  h2  H/faa^)  .  20  log  e 
=  342.9  h2  H/aa 

s 


To  assess  the  bottom  bounce  component  of  noise  from  a  distributed 
surface  source  having  the  above  transmission  characteristics  the 
contribution  from  annuli  are  integrated  (eg  figure  6). 

The  area  dA  of  the  annulus  at  range  r  is 


dA  =  2trr  dr 


Thus  the  total  noise  contributed  by  the  annulus  is 


diDD  =  2irr  .  s  .  k  r  5/2  exp(-a'r) 

DD  O 


Thus  the  integrated  noise  contribution  between  ranges  r j ,  r2  is 


r  2 


J 


BB  =  2  it  so  k  j  exp(-a'r)  .  r  5/ 2  dr 


r  i 


Let 


Then 


2  » 
x  =  a  r 


x3 

f 


iDD  =  4‘IT/a,  s  k  I  exp(-x2)/x2  dx 

DD  O  1 
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where 


Xj  =  /a ' r j 


and  x2 


=  /a  '  r , 


Then 


x2 


1  =  Uv/aJ  s  k  I  l/x2[l-x2  +  xfc/2!  -  x6/3!  +  ...]  dx 

BB  o  j 


*2 


t 


-  4Tf%/a"^  sk  I  (l/x2-l)dx 


if  (x2<lj 


=  4ir/cT^  s  k  [1/Xi  -  l/x2  +  Xj  -  x2 


Hence  taking  logs  and  evaluating  all  constants,  this  becomes: 


IBBfdB)  =  Sq  -  26.67  +  10  log  Ch2  H/cts)  +  5  log  a' 

+  10  log  [  1/i/a  '  r  i  -  l//a'  r2  +  /a  '  r  2  -  /a  '  r  2  ]  (18) 

where  S  =  S  +  60 . 
o 

The  values  of  rt  and  r2  are  determined  by  simple  geometry  and  test 

calculations  as  ri  =  (2H  -  h^j/tan  0^  and  r2  -  100  km.  The  angle  0^  is 

the  critical  angle  defining  the  boundary  between  the  BB  and  direct  path 
components  shown  in  figures  1  and  2.  The  receiver  depth  is  h^Cm).  At 

high  frequencies  x2  becomes  greater  than  1  so  the  approximation  to  the 

integral  is  not  accurate.  However,  tests  show  that  the  monopole 

contribution  given  in  the  next  section  dominates  in  this  case. 

2.7.2  The  monopole  contribution 

Comparison  between  model  predictions  based  solely  on  the  dipoie 

contribution  above  and  noise  measurements  show  major  discrepancies  at 
low  frequency.  The  dipole-only  model  predicts  much  lower  values  than 
those  observed.  In  common  with  other  workers  leg  Wagstaf f ( ref . 25) , 
Wilson! ref . 8) ) ,  this  characteristic  has  been  ascribed  to  the  fact  that 
there  is  not  complete  cancellation  between  direct  and  surface  reflected 
components  contributing  to  the  dipole.  While  there  are  several  ways  of 
rectifying  this,  the  appropriate  approach  still  cannot  be  suggested 
since  actual  physical  mechanisms  are  not  yet  confirmed.  However  the 

interim  solution  adopted  in  the  DUNES  model  is  to  add  a  monopole 

contribution  to  the  dipole  contribution.  This  fills  in  the  dipole  null 

and  implies  that  the  source  mechanism  has  two  components. 
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The  monopole  contribution  is  assumed  to  have  transmission  loss  as  in 
equation  (13)  which  applies  to  propagation  in  the  model  stripping  region 
for  a  deep  source.  This  expression  therefore  describes  transmission 
without  dipole  effects  -  viz 

1  =  —  r  3/ 2  .  exp( -a ' r ) 

m  a 

Integrating  the  intensity  of  an  annulus  of  surface  energy  in  the  same 
manner  as  the  dipole  contribution  yields 


di  =  2irr  .  s  .  r  3/ 2  exp(-a'r)  /  a  dr 

mono  m 


The  integrated  intensity  is 


mono 


r2 


it  s  /a  e  //r  .  dr 


let 


x2  =  a’r,  xi  =  a'ri  and  x\  =  a  r2 


then 


i  =  4tts  /(av^a"*")  I  e  X  dx 


mono  m 


so 


i  =  2tts  /iT/(a/a' ferf(x2)  -  erf(x2)] 

mono  m 

Taking  logs  and  evaluating  constants, 

I  =  10.46  +  S  -  63  -  5  log  a' 

mono  m 

+  10  log  [erf (/r2of )  -  erf (/ r 2a ' ) ]  (19) 

(Note  a'  =  0.23  a  where  a  is  in  decibels/kilometres).  S  is  the 

m 

monopole  source  level  ie  Sq  -  10  dB.  Again  rj  is  given  by  geometry. 
The  value  of  r2  is  taken  to  be  5000  km  unless  there  is  bathymetric 


15 


WSRL-TN-34/89 


blockage.  Estimates  of  the  bottom  bounce  component  involve  the 
evaluation  of  dipole  (equation  (18))  and  monopole  feauation  (19)) 
components.  They  are  combined  by  power  summation.  Under  average 
conditions,  these  are  graphed  in  figure  11. 

2.8  The  extra  down-going  components 

2.8.1  The  bottom  bounce  sector 

Within  the  same  angle  sector  as  bottom  bounce  energy  there  is  an 
additional  component  which  arrives  in  the  downward  direction  only.  This 
component  is  received  directly  from  the  ocean  surface  from  an  annulus 

defined  by  the  intersection  at  the  surface  of  angles  0X  =  cos  1  (C  /Cn ) 

O  B 

and  02  =  cos  1  (C  /C.„)  as  in  figure  2.  For  the  case  of  C  '■  C  then 
O  BB  O  B 

0i  =  o. 

Since  this  component  does  not  have  bottom  interactions  it  is  best 
described  in  the  same  way  as  RSR  paths,  namely  by  equation  112). 

However,  in  this  case  a  more  general  form  of  the  dipole  loss  term  is 
required  to  recognise  a  non-zero  lower  angle  limit  to  the  integration  of 
the  dipole  function.  The  equivalent  expression  to  equation  (5)  for  the 


Ddown  =  “  '  (sin  •+itz02  *  sin  witz©  i )  /  (2itz(02  -  0  j. )  >  (20) 

This  dipole  loss  is  in  addition  to  simple  cylindrical  spreading  from  the 
source  annulus  to  the  receiver.  Using  the  same  expression  for  received 
noise  ( I  <iown  J  as  was  developed  for  the  RSR  case  (cylindrical  spreading) 

we  obtain  (from  equation  112)) 


I 


down 


S,  +8-70-10  log  a’  +  10  log  D, 
k  down 


+  10  log  (exp(-a’rQ)  -  exp(-a’ri)) 


(21) 


The  ranges  of  integration  (r^,  rj)  can  be  estimated  from  the  geometry  of 

the  situation  as  rQ  =  h^/tan  02  and  rt  =  h^/tan  0j.  These  are  typically 

300  m  and  800  m  respectively  for  the  case  of  a  200  m  receiver  depth. 
The  corresponding  value  of  the  downward  component  described  by 
equation  (21)  is  shown  graphically  in  figure  9(a)  as  a  function  of  wind 
speed  and  frequency.  By  comparison,  the  bottom  bounce  component  is 
shown  in  figure  91b).  It  can  be  seen  from  these  figures  that  the 
downward  component  is  insignificant  compared  with  the  bottom  bounce 
component  at  frequencies  below  1  kHz  and  therefore  does  not  affect 
estimates  of  noise  in  the  relevant  vertical  angles.  However,  for  the 
case  of  a  deep  ocean,  deep  receiver  at  higher  frequencies  this  extra 
component  can  be  important  and  it  is  included  in  DUNES.  The  effects  are 
seen  in  the  examples  in  Appendix  I. 
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2.8.2  The  RSR  sector 

There  should  be  an  additional  downgoing  component  included  in  the  RSR 
angle  sector  between  0!  and  0  .  However  this  is  not  currently  included 

in  DUNES  2.3.  Again  it  will  only  be  important  at  high  frequencies. 

2.9  The  direct  path  component 

The  final  wind  generated  noise  component  considered  within  the  DUNES  model 
is  that  arriving  directly  from  the  surface  at  angles  above  cos  1  (C  /C  ) 

O  DD 

or  an  equivalent  critical  angle  (see  figure  2).  An  estimate  which  relates 
to  this  has  been  offered  by  Burgess  and  Kewley (ref . 5 ) .  Their  approach, 
based  on  an  energy  flux  argument,  results  in  a  noise  level  N,  where, 

N  =  S  +  AB  +  8  (22) 

where 


AB  =  10  log  ((b  +  l)/(b  -  1))  (23) 

b  =  bottom  reflection  loss  at  vertical  incidence 

=  antilogjo  (B/10j 

and  3  =  bottom  loss  in  decibels 

This  noise  estimate  assumes  in  essence  that  the  direct  path  contribution 
defined  above  extends  across  all  vertical  angles.  It  is  necessary  to 
remove  from  this  estimate  noise  outside  the  angles  of  interest  -  namely 

noise  within  ±0  =  cos  1  (C  /C  )  .  Since  energy  coming  at  the  receiver 

L  O  DD 

from  angles  below  horizontal  have  one  more  bottom  interaction  than  those 
incident  from  above,  the  energy  N  in  equation  (22)  has  to  be  separately 
apportioned  into  up  and  down  components.  Some  discussion  of  this  is  given 
in  reference  5. 

Based  on  this,  it  can  be  shown  that  the  up  and  down  components  of  direct 
path  noise  are,  respectively 


I  —  S  +  8  +  10  log  a  +10  log  E 
up  °  u  ° 


(24) 


Jdown  =  s  +  8  +  10  loS  ad  +  10  log  E  l25) 

where  E  =  1  -  sin  (cos  1  (C  /C  )} 

O  DD 

a^  =  exp(-za')  exp( - (H-z)2a ' )/ (b-exp( -2Ha ' ) ) 
a^  =  b  exp(-za' )/(b-exp(-2Ha' )) 
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z  =  the  receiver  depth  (km) 
a'  =  the  water  column  attenuation  (nepers/km) 
H  =  water  depth  (km) 

In  the  limits  of  low  frequency  and  no  attenuation 


au  =  1/ (b  -  1) 
ad  =  b/(b  -  1) 


which  are  the  limits  used  to  derive  equation  (23).  In  the  limit  of  b  =  1 
the  above  expressions  allow  a^  and  a^  to  stay  finite.  The  expressions  also 

reduce  the  amplification  for  high  frequencies  by  including  attenuation. 


The  values  of  vertical  incidence  bottom  loss  used  in  the  DUNES  model  are 
derived  from  measurements  taken  from  Burgess  and  Kewley ( ref . 5 ) .  Figure  10 
(based  on  figure  11  in  reference  5)  shows  the  frequency  dependence  of 
measured  vertical  incidence  bottom  loss  at  various  sites  in  the  southern 
hemisphere.  The  trend  is  simulated  by  the  modified  sinusoid  which  is  also 
plotted  on  the  figure. 


The  sinusoid  has  the  form 


6  =  2.9+2  sin  (3.5  log  f  -  6.82)  (dB) 


(26) 


The  sinusoid  is  modified  by  fixing  the  predicted  value  below  32  Hz  at  0.9 
and  above  500  Hz  at  3.88. 

Predicted  levels  of  the  direct  path  wind  noise  component  are  shown  in 
figure  11.  A  correction  factor  is  also  shown  which  accounts  for  the  effect 
of  moving  the  simulated  normal  incidence  bottom  loss  (figure  10)  up  and 
down  to  high  and  low  loss  extremes.  The  final  minimum  loss  is  set  at 
0.5  dB. 

2.10  The  shipping  component 

Although  the  DUNES  model  has  been  motivated  by  a  need  to  predict  and 
understand  wind  generated  noise  it  also  contains  estimates  of  the  shipping 
contribution.  In  most  southern  hemisphere  locations  of  interest,  this 
contribution  is  difficult  to  predict  because  the  ship  densities  are  low  and 
good  historical  information  not  available.  The  shipping  contribution  is 
commonly  the  result  of  a  few  local  vessels  whose  position  cannot  be 
predicted  easily  and  whose  source  spectrum  is  not  known.  There  is, 
nevertheless  a  clear  need  to  include  ships  since  in  some  instances  their 
position  is  known.  There  is  also  interest  in  assessing  the  relative 
contribution  of  ships  compared  to  wind  generated  noise  at  various  locations 
and  frequencies.  For  these  reasons,  shipping  has  been  included  even  though 
in  some  instances  with  less  detailed  description  than  corresponding  wind 
noise  components. 
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2.10.1  Source  spectrum 

One  of  the  greatest  problems  in  assessing  the  ship  component  is  not 
knowing  an  appropriate  source  spectrum.  The  source  level  of  a  ship 
depends  upon  ship  speed,  horsepower,  mechanical  condition  and  the  sea 
state.  These  factors  vary  in  a  way  which  is  hard  to  predict.  Under 
conditions  of  high  ship  density,  some  advantage  can  be  taken  of 
generalised  ship  statistics  where  individual  ship  variation  is  absorbed 
in  the  overall  picture.  This  approach  has  been  adopted  by  some 
authors ( ref . 26 , 27 )  and  has  greatest  relevance  in  high  shipping  density 
areas.  In  regions  where  ship  noise  is  controlled  by  a  few  individual 
ships,  a  more  specific  prediction  based  on  noise  from  ship  classes  is 
more  desirable.  In  most  instances,  however,  specific  ship  passages 
cannot  be  predicted  with  sufficient  precision  and  probabilistic 
solutions  must  be  sought.  To  this  end  some  initial  surveys  were  taken 
of  ship  traffic  around  New  Zea land( ref . 28)  .  The  general  distribution  of 
ship  length  (which  correlates  broadly  with  source  level)  was  shown  to  be 
similar  to  world  wide  distributions  (figure  12).  The  most  likely  ship 
length  is  around  600  ft  which  applies  to  general  cargo  vessels. 
Heine ( ref . 27 )  has  reported  on  measurements  of  such  a  vessel  and  these 
are  summarised  in  figure  13  for  the  broadband  radiation.  Specific 
frequency  lines  will  bo  up  to  10  dB  higher  than  the  continuous  spectrum 
but  these  are  not  the  subject  of  the  noise  prediction.  Line  clutter  is 
regarded  as  easily  distinguishable  from  background  noise  and  falls 
within  the  field  of  target  identification. 

The  general  shape  of  the  continuous  spectrum  (figure  1 3  ('  b  j  )  car.  be 
described  by  three  parameters  -  low  frequency  slope,  high  frequency 
slope  and  peak  level.  Low  frequency  roll-off  is  determined  by  surrace 
dipole  effects  and  should  be  included  as  part  of  the  transmission  loss 
not  source  level.  However,  this  involves  estimating  effective  source 
depth  of  ship  radiation  and  a  level  of  investigation  beyond  the  needs  of 
the  DUNES  model.  The  low  frequency  roll-off  is  retained  as  6  dB/octave 
in  the  model.  The  high  frequency  slope  (15  dB/ octave  i  has  been  found  to 
be  too  steep  and  ^o  10  dB/octave  is  in  the  model.  Measured  shipping 
noise  examples  are  seen  in  reference  29  (figures  2  to  6,  2  to  9  and 

2  to  14).  Maximum  source  level  is  an  input  parameter  in  DUNES.  The 
default  value  of  167  dB  (as  in  figure  13(b))  seems  to  produce  acceptable 
noise  predictions  under  conditions  where  ship  noise  predominates. 

2.10.2  Enhanced  ships 

It  is  well  established  that  sources  over  distant  bathymetric  slopes  or 
at  high  latitudes  give  rise  to  enhanced  noise  contributions.  This  was 
discussed  in  connection  with  wind  generated  noise  in  Section  2.3. 

Ship  sources  in  such  situations  are  dealt  with  in  either  of  two  ways 
depending  whether  it  is  easier  to  describe  them  as  individual  vessels  or 
as  a  continuous  ship  lane. 

The  noise  contribution  from  individual  ships  ( L,,,. )  is  estimated  using  a 

orl 

cylindrical  spreading  transmission  loss,  viz 


L  „  =  S  -  70  -  10  log  R  -  aR 

oh 


(  2"  ) 


where  S 
R 

a 


ship  source  leve  1  at  a  given  frequency- 
range  to  the  receiver  (km) 
water  column  attenuation  (dB/km) 
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If  the  shipping  distribution  is  best  described  as  a  continuous  lane, 
then  a  modified  form  of  equation  (3)  (Section  2.3)  is  used. 

The  noise  contribution  (Icr)  of  a  ship  lane  at  a  distance  R  from  a 

b  J_i 

receiver  is 

I$L  =  SA  -  70  -  3  +  K  -  a  R  /xi/4R2  +  1  (28) 

where  S. 

A 

S 

S\ 

Xi 

R 
K 

a 

The  noise  contribution  predicted  by  equations  (27)  or  (28)  come  by  way 
of  RR  paths  as  defined  in  figure  2. 

2.10.3  Local  ships 

The  transmission  of  noise  from  local  ships  is  via  RSR.  bottom  bounce  and 
direct  paths  in  a  way  analogous  to  wind  generated  noise  (Sections  2.-+  to 
2.91.  This  approach,  however,  involves  the  complication  of  needing  to 
know  the  source  depth  of  ship  noise  radiation  which  is  difficult  to 
estimate.  Consistent  with  the  approach  described  in  Section  2.10.2,  a 
single  transmission  loss  is  used,  wnich  relates  only  to  bottom  bounce 
modes.  It  is  found  that  spherical  spreading  transmission  loss  plus 
attenuation  provides  a  reasonable  approximation  over  the  shipping 
frequency  band  for  shallow  sources.  While  this  description  does  not 
allow  for  different  bottom  loss  values  it  represents  the  current 
implementation  in  DUNES  and  has  been  found  to  give  useful  results. 

Accordingly,  the  case  of  discrete  ships  is  handled  as  follows.  The 

received  level  (I_„)  is  given  by  receiver  is 
bn 

lc„  =  S  -  60  -  20  log  R  -  cxR  (29) 

bn 

whore  S  =  ship  source  level  (Section  2.10.1) 

R  =  range  to  receiver  (km) 
a  =  water  column  attenuation  (dB/km) 

When  local  shipping  is  most  easily  descrxDed  as  a  uniform  shipping  lane, 
then  the  noise  contribution  is  approximated  by 


=  source  level  of  shipping  per  kilometres 
=  S  +  10  log  (SN/X!) 

=  single  ship  source  level  (dB)  re  1  yPa2/Hz  at  1  km) 
=  number  of  ships  on  lane 

=  half  the  length  of  the  symmetrical  ship  lane  (km) 

=  CPA  range  of  lane  to  receiver  (km) 

=  10  log  [InfXi/R  +  /xf/R2  +  1)] 

=  water  column  attenuation 
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where  parameters  are  the  same  as  for  equation  (28).  The  local  ship 
noise  contribution  calculated  from  either  equations  (29)  or  (30)  is 
assumed  to  cover  the  vertical  angle  range  relevant  to  bottom  bounce 
transmission  as  defined  in  figure  2.  There  is  no  calculation  of  the  RSR 
(convergence  zone)  contribution  since  this  requires  information  which  is 
too  specific  geographically.  The  level  of  the  bottom  bounce 
contribution  is  extended  to  cover  the  RSR  vertical  angle  range 
(figure  2)  to  avoid  discontinuities.  If  a  ship  is  known  to  be  present 
and  within  a  convergence  zone,  its  contribution  can  be  estimated  using 
equation  (27).  This  is  relevant  since  transmission  loss  from  a 
convergence  zone  follows  approximately  a  cylindrical  spreading  law. 

2.11  Storm  noise  component 

Following  the  development  of  shipping  lane  contributions  to  the  noise,  a 
component  of  wind  generated  noise  due  to  non  local  storms  is  easily- 
formulated.  Long  distance  high  latitude  storms  and  coastal  enhanced  storms 
are  readily  accounted  in  the  previous  sections  while  a  local  storm  is 
incorporated  in  the  direct  path  component.  The  storm  is  modelled  as  a 
finite  width  lane  transmitted  via  bottom  bounce  modes  in  the  same  manner  as 
local  ships  in  Section  2.10.3. 

Thus  the  received  level  due  to  storms  is  given  by 


I  =  S  +  10  log  V  -  60  +  3  -  10  log  R  +  10  log  (tan'^/R) 

o  1 

-aR/l  +  Xj/iR2 


where  S 
R 
X! 


source  level  of  storm  noise  per  kilometres2 

range  to  storm  (km) 

half  the  length  of  storm  (km) 


W  =  width  of  storm  (km) 


Follow  the  local  ship  procedure,  the  noise  contribution  is  assumed  to  cover 
the  vertical  angle  range  relevant  to  bottom  bounce  transmission  and 
extended  to  the  RSR  range  to  avoid  discontinuities.  Locating  a  storm  over 
a  known  convergence  zone  region  can  be  simulated  by  using  the  slope 
enhanced  description  as  the  propagation  law  is  then  more  appropriate 
however  the  angular  spread  will  be  wrong. 

2.12  Receiver  depth  effects 


It  has  been  observed  that  the  location  of  the  receiver  close  to  the  sea 
surface  can  cause  a  reduction  in  noise  level  due  to  surface  decoupling 
effects.  Thus  the  components  of  noise  via  BB,  RSR,  and  RR  paths  will  be 
increasingly  sensitive  to  the  depth  of  the  receiver.  The  following 
reduction  factors  are  used  to  account  for  this  effect. 


First  the  decoupling  depth  (in  wavelengths)  for  each  component  is  found 
us  ing 
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dRR  »  0.25/si„(0RR/2) 

dRSR  =  °-2:>/sln^o-5(^0RR  +  0Rsr}) 

dBB  =  °-25/sin(0-5(6BB  +  max(°RR’°RSR))) 


where 

il 

OS 

c* 

© 

cos  1 

(Co/cs) 

SRSR 

cos  1 

(W 

0BB 

cos  1 

(Co/CBB> 

These 

depths 

correspond  to  the  first  maximum  for  the  Lloyd  mirror 

interference  pattern.  The  average  arrival  angle  for  each  noise  component 
is  used  to  approximate  the  typical  decoupling  depth. 

If  the  receiver  depth  (in  wavelengths)  is  less  than  these  decoupling  depths 
then  the  reduction  factors  are  given  by: 

fRR  =  sir>2(2lTZ/(4dRR)) 

fRSR  =  sln2(2lTZ/(4dRsR)) 
fB0  =  sin2  (2ttz/  (4dgB) ) 


where  z  is  the  receiver  depth  in  wave  lengths ( ref . 30 )  .  More  precise 
estimates  can  be  made  using  the  actual  sound  speed  prof i les ( ref . 30 )  rather 
than  the  first  order  isospeed  assumptions  used  here. 


3.  RUNNING  THE  FORTRAN  PROGRAM 


3 . 1  General 

In  its  present  form  the  model  is  set  up  to  operate  either  interactively 
from  a  terminal  or  from  a  file  of  input  parameters.  In  either  case,  the 
main  issue  is  an  appropriate  definition  of  input  values  and  this  is  the 
point  to  be  discussed  in  this  part  of  the  report. 

The  inputs  are  generally  not  required  to  be  specified  with  great  precision. 
The  concept  of  the  model  is  in  terms  of  generalities  and  averages  so  the 
inputs  can  reflect  this  philosophy.  The  decibel  law  allows  for 
considerable  input  latitude  before  output  levels  (dB)  are  affected  greatly. 

Appendix  I  reproduces  a  typical  terminal  session  annotated  with  additional 
components.  In  what  follows  in  this  part  of  the  report,  the  individual 
input  parameters  and  options  are  discussed  with  the  aim  of  assisting  a  user 
to  define  input  values  appropriate  to  a  specific  modelling  requirement. 

3.2  Parameter  inputs 

Reference  to  Appendix  I  will  show  that  the  program  inputs  can  be 
conveniently  grouped  into:  initial  parameters,  environmental,  local  noise, 
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distant  noise  and  shipping  inputs  graph  format  options  and  output  data 
options.  These  separate  sets  of  inputs  are  discussed  in  following 
Sect  ions . 

It  has  been  found  convenient  to  use  an  input  worksheet  to  assist  in 
organising  parameters  relevant  to  a  given  site  and  situation.  The 
worksheet  is  shown  in  Appendix  II  in  blank  form  and  in  completed  form  with 
data  used  in  the  test  run  to  be  described. 

3.2.1  Initial  parameters 

The  implementation  at  WSRL  on  an  IBM  3033  is  initiated  by  a  CLIST 
command  which  is  set  up  to  run  the  program.  For  the  VAX  and  IBM  PC 
versions  similar  DCL  and  BAT  files  can  be  used,  respectively. 

3.2.2  Environmental  inputs 

The  next  set  of  inputs  (Appendix  I)  define  the  environment. 

(a)  Ocean  area 

This  is  used  to  select  from  data  within  the  program  appropriate 
parameters  for  the  high  latitude  wind  contribution  and  sound 
attenuation.  These  are  found  in  subroutines  HLAT  and  ATTENU. 

(b)  Receiver  latitude 

Used  in  subroutine  HLAT  to  determine  distance  from  high  latitude 
winds.  It  is  also  used  to  distinguish  the  North  and  South  Atlantic 
and  Pacific  Oceans  in  ATTENU. 

U  c )  Season 

Also  used  only  in  HLAT  in  initial  parameter  definition.  The  two 
seasons  recognised  in  the  program  are  southern  summer  (February)  and 
winter  (August ) . 

(d)  Receiver  depth  (metres) 

This  is  used  in  subroutine  LOCAL  in  the  area  used  to  estimate  the 
bottom  bounce  component.  The  specific  significance  of  receiver  depth 
is  in  the  calculation  of  the  initial  range  of  integration  (rj  in 

equation  (18)).  This  is  the  range  from  the  receiver  which 
corresponds  to  the  surface  interaction  of  the  ray  which  just  grazes 
the  bottom  and  passes  through  the  receiver.  The  depth  is  also  used 
for  decoupling  factors. 

(e)  Sound  speeds  (m/sl  and  water  depth  (km) 

In  order  to  calculate  the  vertical  angle  limits  which  correspond  to 
various  noise  components,  sound  speeds  at  the  receiver  location  need 
to  be  specified  (see  figure  2).  Parameters  required  are  sound  speeds 
at  the  surface,  receiver  depth  and  bottom.  The  water  depth  also 
needs  to  be  specified. 

(f)  Range  of  transmission  blockage 

At  low  frequency,  the  RSR  component  can  be  affected  by  surface  noise 
to  many  thousands  of  kilometres  (figure  7).  This  distant  component 
would,  however,  be  removed  if  bathymetric  or  land  mass  blockage 
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intervened.  To  allow  for  this  possibility,  the  range  to  significant 
transmission  blockage  (in  thousands  of  kilometres)  is  specified  in 
each  of  four  quadrants  if  near  to  a  coastal  margin,  for  example 
2  quadrants  could  be  set  to  0.1  thousand  kilometres  (100  km)  the 
remaining  2  quadrants  to  the  default  value  of  5  (5000  km)  to  indicate 
no  blockage.  The  quadrants  are  centred  upon  the  NSEV  directions  to 
provide  azimuthal  variation. 

3.2.3  Local  wind  noise  inputs 

(a)  Slope  of  bottom  loss  versus  angle 

Bottom  loss  is  used  in  the  model  as  the  slope  of  the  bottom  loss 
versus  angle  (equation  (14)).  This  is  calculated  internally  as  a 
function  of  frequency  in  a  way  which  seems  to  give  good  agreement 
with  average  transmission  loss  estimates.  The  mode)  input  is  a 
scalar  which  is  a  simple  multiplier  of  the  decibels  per  radian 
average  bottom  loss  estimate.  A  multiplier  greater  than  one  applies 
to  high  loss  areas  and  less  than  one  to  low  loss  areas. 

(b)  Vertical  bottom  loss 

The  vertical  bottom  loss  is  used  to  assess  direct  path  local  wind 
noise  (equations  (22),  (23)).  The  internal  calculation  is  based  on 
figure  10.  The  input  control  over  vertical  bottom  loss  is  simply  to 
move  the  average  curve  (figure  10)  up  or  down  a  specified  number  of 
decibels.  The  range  of  measured  values  shown  on  figure  10  implies 
that  input  limits  of  ±1.2  dB  covers  the  range.  It  should  be  noted 
that  the  model  limits  internally  the  vertical  bottom  loss  to  be 
greater  than  0.5  dB  for  any  input  conditions. 

(c)  Local  noise  angle  limits 

Energy  received  within  vertical  angles  which  correspond  to  rays 
hitting  the  sea  floor  above  the  critical  angle  is  treated  separately 
from  the  bottom  bounce  component  and  is  termed  'local'  noise.  The 
angle  limits  which  apply  to  this  component  are  between  the  vertical 
(up  and  down)  and  an  angle  specified  in  either  of  three  ways  as 
program  input.  Firstly,  the  angle  can  be  specified  directly. 
Suggested  values  are  the  angle  with  bottom  loss  -3  dB  down  from  the 
90°  value  or  the  angle  of  intromission  if  C  /C  >  1.  Secondlv,  it 

ODD 

can  be  calculated  from  the  bottom  critical  angle  if  the  sub-bottom 
sound  speed  is  given.  Thirdly,  a  default  condition  assumes  the 
sub-bottom  sound  speed  is  10%  greater  than  the  water  sound  speed 
(previous  input)  at  the  sea  floor. 

(d  >  Help  sheet 

Appendix  III  shows  a  help  sheet  to  use  when  determining  what  angles 
and  losses  to  use  for  the  local  wind  noise  inputs. 

3.2.4  Distant  wind  noise  inputs 


The  distant  wind  noise  components  are  defined  here  as  originating  at 
high  latitudes  on  bathymetric  slopes.  These  components  are  assumed  to 
be  RR  transmission  types  (Section  2.4). 
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(a)  Fraction  of  high  latitude  wind  and  their  bearings 

This  input  allows  the  user  to  enter  an  estimate  of  the  fraction  of 
the  high  latitude  wind  band  which  is  acoustically  visible  from  a 
given  receiver  location.  This  is  intended  to  provide  the  means  of 
dealing  with  bathymetric  or  land  mass  blockage.  The  bearings  are 
usually  0°  and  180°  but  can  be  changed  to  allow  for  blockage  in  those 
directions . 

(b)  Slope  enhanced  w’ind  lanes  or  storms 

One  option  here  is  tu  specify  no  wind  lanes,  however  a  default  lane 
is  available  as  a  rough  average  of  the  contribution  from  wind  noise 
enhanced  by  bathymetric  slopes.  The  specification  of  the  default 
lane  is  1500  km  long,  100  km  wide  at  a  distance  of  500  km  and  270° 
from  the  receiver.  A  12  kn  wind  is  assumed  to  be  blowing  over  the 
region.  Alternatively,  wind  lane  or  storm  data  above  can  be  entered 
for  up  to  5  independent  slope  enhanced  or  storm  regions.  Their 
bearings  provide  the  azimuthal  variation.  Note  that  the  lanes  are 
assumed  to  be  symmetrical  about  their  point  of  closest  approach. 

3.2.5  Shipping  inputs 

Up  to  five  independent  shipping  inputs  may  be  specified  (or  no  inputs  if 
required) .  Shipping  inputs  are  in  two  forms  -  either  as  a  shipping  lane 
or  as  individual  discrete  ships.  In  each  of  these  categories  the  ships 
can  be  specified  as  enhanced  or  not.  The  enhanced  option  implies 
transmission  loss  relevant  to  ideal  down  slope  propagation  over 
bathymetric  features.  This  option  (for  discrete  ships /  is  also  suitable 
for  a  vessel  known  to  be  in  a  convergence  zone. 

If  a  ship  iane  is  specified  then  lane  length  and  total  number  of  vessels 
on  the  lane  is  required.  Note  that  the  lane  is  assumed  to  be  symmetrical 
about  its  point  of  closest  approach.  In  all  cases  of  ship  input  the 
source  spectrum  is  specified  by  a  peak  level  and  a  frequency  at  which 
the  upper  roll-off  (10  dB/octave)  begins.  The  bearings  are  needed  for 
azimuthal  variation. 

3.2.6  Output  graphics  options 

There  are  two  output  wind  speed  options  -  either  to  specify  a  wind  speed 
(5  to  40  kn)  for  full  plot  of  all  noise  components  at  that  wind  speed  - 
or  to  request  total  omnidirectional  noise  at  all  wind  speeds  to  be 
shown.  These  two  options  are  shown  in  Appendix  I  ig  and  i).  There  are 
also  plots  of  vertical  and  horizontal  noise  directionality  available  on 
request . 

The  units  of  the  vertical 
omnidirectional  energy  fdB 
(dB  re  1  yPa2/Hz  -  srf  The 
models  to  compare  with  data. 

3.2.7  Output  data  options 

Tables  of  omnidirectional,  vertical,  horizontal  and  three  dimensional 
noise  are  available  to  be  output  before  DUNES  finishes.  The  first  three 
files  are  Formatted  while  the  last  is  Unformatted.  The  outputs  are  in 
columns  of  two  as  in  the  following  examples 


noise  can  be  specified  as  either  total 
re  1  yPa2/Hzj,  or  energy  per  steradian 
latter  units  are  used  in  most  prediction 
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N  Angle - 
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leve  1 

3.3  Final  notes 

Some  notes  on  IBM  PC  usage  are  provided  in  Appendix  IV.  If  errors  or 
deficiencies  are  found  in  the  DUXES  model  the  authors  would  be  pleased  to 
know  about  them  to  ensure  upgrading  of  the  model  can  continue. 

4.  CONCLUSIONS 

At  present  the  ability  of  the  current  model  to  be  extended  to  lower 
frequencies  and  to  shallow  water  (<200  m)  is  being  investigated.  The 
FORTRAN  77  code  is  available  from  the  authors. 

The  development  of  the  DUNES  model  presented  here  explains  the  methods  used 
and  why  it  has  proven  to  be  an  effective  model  in  giving  realistic  predictions 
of  noise  level  and  directionality.  However,  there  are  a  number  of  assumptions 
made  by  the  authors  regarding  source  depths  and  directionality  and  only 
further  studies  in  these  fields  will  resolve  their  true  forms.  We  believe 
that  DUNES  can  be  used  to  further  our  understanding  of  these  phenomena. 
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APPENDIX  I 

TERMINAL  OUTPUTS  DURING  DUNES  RUN  ON  IBM  AT  VSRL 
1.1  Initial  parameters 


DUNESgS  PLOT  (TEC) 

FILE  FTGSFDG1  IMS  ALLOCATED  AND  HAD  KOI  FREED  DCF ORE  THE  CURRENT  ALLOCATION. 
FILE  FTRSFM1  HAS  KEN  ALLOCATED  SHR. 

FILE  FTMF001  HAS  KEN  ALLOCATED  SHR. 

FILE  FTMFM1  HAS  DEEM  ALLOCATED  SHR. 

FILE  FT10FM1  HAS  SEEN  ALLOCATED  SHR. 

FILE  FTUFM1  HAS  KEH  ALLOCATED  SHR. 


1.2  Environmental  inputs 

DUNES  (2.3)  NOISE  PREDICTION  AO DEL 


GENERAL  ENWIRONNENTAL  INPUTS .  )  Needed  for 

ENTER  OCEAN  AREA(  1  •PACIFIC,  2-IHDIAN,  3-ATLANTIC  )  high  latitude 

J  )  wind  components 

ENTER  RECEIVER  LATITUDE' DECS)  -SOUTH  NEGATIVE 
7 

-30 

ENTER  SEASON  ( 1‘FEMUARV,  2 -AUGUST) 

? 

I 

ENTER  RECEIVER  DEPTH(NJ. .EG  209 
7 

2M 


DEFAULT  SOUNDSPEEDS  (H/S)  ARE) 

1S2S(  SURFACE ).  I  49*<  RECE IVER ),  IS3K  DOTTON ) 

ENTER  1  TO  ALTER  THESE  VALUES 
7 
1 


ENTER  SOUNDSPEED  VALUES  AT . 

SURFACE, RECEIVER. DOTTON  < 1S2D, 1499, 153*) 

IS 27  149*  1539 

SPECIFY  WATER  DEPTH  (KN).. 

7 

4.5 

ENTER  RANGE  (KHS1D9G)  CLEAR  OF  TRANSMISSION  > LOCKAGE  )  Important  for 

IN  4  QUA9RANTS(R,1M,M  AND  27G  ♦-  45  KGS  )  Strong  RSR 

E.G.  S.S.S.5  IN  AID-OCEAN)  )  component 

.1  .1  5  S 
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1.3  Local  wind  noise  inputs 


IWUTS  FOR  local  UINOHOISE  COMPONENT . 

ENTER  SCALE  FACTO*  FOR  SLOPE  Of  BOTTOM  LOSS  VS  ANGLE 

' 1-AVERACE  * 

•> 

l 

ENTER  OB  'VARIATION  FROM  AVERAGE  VERTICAL  BOTTOM  LOSS 
(♦1.2  AND  -1.2  ARE  HIGH/LOU  EXPERIMENTAL  UNITS,  B-AV) 

7 

0 

. TO  KFIHE  LOCAL  NOISE  LINITSi 

ENTER  EITHER  -  ANCLE  LIMIT  FOR  LOCAL  NOISE  (DEO) 

OR  -  SUB-BOTTOM  SOUNDSPEED  <M/S) 

OR  -  ZERO  FOR  DEFAULT 

• 


)  High  loss  >  1 
)  Low  loss  <  1 


)  Default  is  10°o  jump 
)  in  sound  speed  from 
)  water  to  bottom 


1-4  Distant  wind  noise  inputs 


INPUTS  FOR  DISTANT  UINONOISE  COMPONENTS 


ENTER  FRACTION  OF  HIGH  LATITUDE  WIND  REGION 

(NORTH, SOUTH) . E.C.  1,1 

7 

t.s  i 

NOV  ENTER  THE  N  AND  S  SEARINGS,  USUALLY  D,1M  DEG 
0 

•  iso 


ENTER  NUMBER  OF  SLOPE  OR  STORM  WIND LANES 
t-l*ONE  DEFAULT  SLOPE  ENHANCED  LANE;  ••NONE  S-MAX) 

*9 

i 


FOR  UIHDLANE  NUMBER  1  ENTER  FOLLOWING 

LAME  LENGTH  (KM) 

? 

!••• 

IS  THIS  A  SLOPE  ENHANCES  LAM  (V/MIT 
V 

lane  width  <kh> 

■» 

26 

CPA  RANGE  FROM  RECEIVER  (KM) 


BEARING  OP  LAME  FROM  THE  RECEIVER  (KG) 

T 

27  0 

UINSSREED  AT  UINBLANE  (KNOTS) 

•» 

IS 


)  To  cover  land  mass 
)  shadowing  of  high 
)  latitude  winds 


)  Default  lane  is  1500  km 
)  1500  km  long,  100  km  wide, 

)  500  km  range,  at  270°, 

)  enhanced  with  12  kn  winds 
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1.5  Shipping  inputs 


INPUTS  FOP  SHIP  CONTRIBUTION . 

ENTER  NUNBER  OF  SHIP  INPUTS  TO  FOLLOW  <•  TO  S) 

7 

2 


iitimtiiiitittittiii 
FOR  SHIP  INPUT  NUNBER  1  ENTER  FOLLOWING* 

LANE  LENGTH  <KH) 

(LANE  LENGTH  INDICATES  DISCRETE  SHIP  OPTION) 

7 

B 

IS  THIS  SHIP  CONTRIBUTION  SLOPE  ENHANCED?  (V/N>» 

V 

RANGE  Of  SHIP(S)  IN  KR 
7 

1M 

BEARING  OF  SHIP<S>  FRON  RECEIUEP  (DEG) 

? 

71 

NUNBER  OF  SHIPS  AT  RANGE  AND  BEARING 
7 
1 

PEAK  SOURCE  LEUEL  OF  AUER AGE  SHIP  (E.G.16S  DB) 

7 

IRS 

AND  FREQ.  BREAKPOINT  OF  SOURCE  SPECTRUM  (E.G.M  HZ) 


iiiimiiiiittmtmt 
FOR  SHIP  INPUT  HUNGER  2  ENTER  FOLLOWING! 

LANE  LENGTH  <KN) 

‘LANE  LENGTH  *0  INDICATES  DISCRETE  SHIP  OPTION) 

7 

1000 

IS  THIS  SHIP  CONTRIBUTION  SLOPE  ENHANCE D7  <V/M)i 
H 

CPA  RANGE  FRON  SHIP  OR  LANE  TO  RECEIVER  (KN) 


BEARING  OF  SHIP  OR  LANE  FRON  RECEIVER  (DEG) 

7 

270 

AVERAGE  NUNBER  Of  SHIPS  (IF  LANE  TOTAL  NUNBER) 

? 

20 

PEAK  SOURCE  LEVEL  OF  AVERAGE  SHIP  (E.G.I6S  DB) 

7 

1SS 

AND  FREQ.  BREAKPOINT  OF  SOURCE  SPECTRUM  (E.G.M  HZ) 
7 

SB 


Discrete  ship  option 


Ship  lane  option 
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1.6  Output  options 


ENTER  UINDSPEED  FOR  ONNI  PLOT  'KNOTS) 

EITHER  ONE  JAIL1E  BETWEEN  S - 40  OR  6  FOR  FULL  RANGE 

t 

IS 


OUNES  PREDICTION  IN  PROGRESS 


•S/«S/88< 

INPUT  PLOT  SIZE  AN*  DEVICE  UNIT 

■> 

1.1 


)  Dummy  values  tor  IBM;  VAX 
)  and  PC  values  are  not 


1.7  Omnidirectional  noise  at  specified  (kn)  wind  speed  with  all  components 
shown 


DUNES  2.3  OMNIDIRECTI0NRL  NOISE 

MIND  LANE  PARAMETERS 


PACIFIC  OCEAN  -  FEBRUAAT  SHIP  LANE  PARAMETERS 

ninospkc'is.’o  Snots  type  len  pang  shps  lEvl  pern  bear  type  len  mioth  rang  bean  mspo 

BBT. LOSS  GRROlENTi  1.0  SHIP*  0  100  1  185  60  73  SLPE  1000  25  100  270  15 

VERTICAL  BOT.LOSSi  0.0  LANE  1000  200  20  165  60  270 

NINOLRNC  FRACTIONS! 

0.5  IN)  1.0  (S) 

RO  1NH200.0EPTMIH.50  KM 

C3.  CB.  CB:  1327.  1HB6.  1338.  QUAD.  RANGE  N.  3.  E.  N«  100.  100.  3000.  3000.  KN 
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1.8  Vertical  directionality 


CHIT*  UtrriCftl  HOIK  UNITS  (ONNIDIKCTIOHAL/STCMOIAMS)  ON  (XJIT  UCrriCM  HOIK  SECTION. . CO/S'Ol 
NOTE  THAT  FON  UCNTIC6L  MTANMC  OCITFUT  UK  S.  8CF0K  0 
5 

UIN96NCE9  •  15. • 


NOUHOSNCEDS-SUN.RCC.HlN.SOTl  1527.68  1498.68 

AHQLC  LIMITS -SO, SO, L0C>  11.2  13.3 

NN.KSN  SWITCH  FL466  1.868  t.NM 

UNITS  •  0*.  'UM8**/<ie-STt) 


1478.66  1539.66 

27.8 

)  Surface  and  bottom  grazing 
)  angles  and  local  angle  limit 


FNCO 

HU 

l»D 

RSN 

NR 

HIGH) 

H  I  OHO 

OMNI 

28. 

63.5 

63.5 

74.99 

76. *8 

57.57 

58.48 

77.22 

28. 

63.9 

64.6 

75.68 

71.68 

57.47 

58. 38 

77.88 

38. 

64.5 

64.5 

76.12 

72.76 

57.42 

58.33 

78.63 

48. 

65.4 

66.4 

75.48 

73.93 

56.43 

57.54 

79.56 

58. 

65.6 

65.7 

75.57 

74.11 

54.98 

S8.48 

79.65 

84. 

65.2 

66.2 

75.47 

73.88 

52.99 

54.87 

79.24 

88. 

63.3 

63.4 

74.44 

72.65 

58.55 

63.33 

77.71 

188. 

61.6 

61.7 

73.49 

76.35 

48.93 

58.34 

76.18 

188. 

59.9 

68.6 

72.37 

68.88 

47.53 

61.58 

74.43 

188. 

57.8 

57.9 

78.58 

64.88 

46.66 

68.88 

71.63 

888. 

56.4 

58.5 

68. « 

61.16 

48.34 

58.88 

69.31 

888. 

58.3 

55.5 

67.73 

57.48 

46.81 

58.16 

67.39 

318. 

53.8 

54.1 

65.96 

54.13 

45.46 

58.25 

68.63 

488. 

S2.8 

53.1 

64.58 

68.48 

45.88 

58.48 

64.48 

884. 

51.8 

52.2 

83.28 

58.78 

48.82 

58.66 

63.51 

838. 

58.4 

51.6 

61.61 

48.88 

48.48 

58.68 

62.28 

m, 

49.4 

56.1 

88.47 

48.46 

48.35 

58.58 

61.42 

1888. 

43.7 

45.1 

54.43 

38.81 

43.38 

47.89 

56.88 

36.8 

39.1 

45. 98 

3.63 

38.88 

44.15 

51.69 

8488. 

29.8 

34.5 

36.77 

8.15 

31.37 

39.87 

46.44 

01/11/88  DUNES  12.3) 

Msreco  =  is.o  knots  veaticpl  noise  LEVELS  D0//UPfl»»2/ (SR-HZ) 


=10  08  SCALE 

■  DOWN 


46  _L 


98  -  7  J 

u 

L _ / 

94. 4  i 

[ 

E 

n — > - 

L _ 

-90  -70 


-SO  -90  -10  10  30 

BEAN  ELEVATION  I  DEC) 


HE. 

HZ 

6400 

OMNI 

46.4 

44.2 

3200 

51.7 

47.8 

1600 

56.9 

SO.  6 

800 

61.4 

SO. 6 

038 

62.3 

50.7 

304 

63. S 

SO. 4 

400 

64.5 

50.2 

316 

65.5 

SO. 2 

252 

67.4 

.50.2 

200 

69.3 

.50.6 

160 

71.6 

.51. 5 

126 

74.4 

11:1 

18° 

7?V 
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1.9  Horizontal  directionality 


rNTFP  -£QTICAL  NOISE  UNITS  (OHNlPISECTIONAl.'STERADlAHS  >  OR  011  IT  OEPTICAL  HOISt  SECTION. 
NC'E  T-iT  FOP  "ERTICAl  DATAIaJE  OltPUT  LHC  5.  IEF0RE  0 

"to  *aNT  HORIZONTAL  NOISE  ROSES  >  Y/N  )R 

UH->*  fPEOuENCIES  FOR  PLOTTING  .  10-10KM2  ’ 

Input  6  OAUJCS,  NEAREST  1-3  OOTAuES  USEP,  P  OK 

•> 

50  508  POP# 

FREQUENCY  'HZ',  REOUESTED  (HZ' 

sp.eoooeoo  w.iinm 

UIND  SPEED  (KN ) 

15.0WMM 

HORIZONTAL  NOISE  SOURCES  M//UPA*«2:;<HZ-DCS> 

I*S  SECTOR  0EAPIHG(DCG '  AM  UIDTMlDEG) 

ncrt*-  and  south  hilat  noise 


4?.?'94l59 
s’.ei36938 


P.__ 

180. 


-9-i 


NORTH,  SOUTH,  EAST,  UEST  ♦  8ASIN 

43.SIHI36  - 

43. S 114136 
S0. 0-76367 

50.0-76367 
UIHD  LANES 
41,7400818 
SHIP  LANES 
74.9943848  73, 

50, 1631927  870 


0.1 - 

180.000000 

90.0000000 

270.000000 


270. 


9.00700062 

55.2399832 

Qi  iiMMi 

90. 

90.1 
90.1 


78.7100034 


08.8442932 


LIST  OF  HORIZONTAL  NOISE  LEVELS  AND  THEIR  AMLES 
DECREES,  08/AlPA»«a/(HZ-0E0) 


1  « 


5O.862KS0 
43.5114130 
50.077630T 
75.0003000 
50.0776387 
43.5114130 

57. 4046921 

57.4046921 
43.5114130 
59.9776397 
59.6719906 
53.4363687 
50.6719666 

59.9770307 
43.5114136 

..  50.8626556 

OMNIDIRECTIONAL  NOISE  D8^/UPA**2/HZ 
■>9.6499176  15.9096690  KNOTS 


2 

4.53364931 

3 

45.I900900 

4 

72.5000000 

5 

73.5000000 

6 

135.000000 

7 

152.380508 

8 

180.000000 

9 

207.619498 

10 

225.999999 

11 

230.641998 

12 

235.877983 

13 

304.122079 

14 

309.367910 

15 

315.000000 

16 

355.466300 

01/11/88  NOISE  ROSE 


DUNES  (2.31  HI  NO  SPEED  =  15.0  XN  LEVELS  0B//UPIUi2/ (HZ-DEC) 

0, 


FREQUENCY  (HZ) 
50.0 

OMNI  NOISE 

DB//UPfl»»2/HZ 

79.6 


(0'5'0> 
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1.9  Horizontal  directionality(Contd. ) . 


FRCOUEHCV  (HZ),  WOOCS TO  (HZ) 


UINC  Mn  (KM) 

ll.MNMt 

HORIZONTAL  HO  I  SC  tOURCti  D»//UPAUa/(HZ-DCS ) 

ITS  SCCTOA  lEARIMCCDCG)  AND  UtDTH(DCC) 

NORTH  AND  SOUTH  HILAT  NOISE 
-9. $7468128  9.96799262 

-19.2733612  180.900009  55.  2319932 

north,  south,  east,  UCST  ♦  IASIN 
35.9896627  9. 999999999C*99  99.9999999 

IS. 988662?  199.999999  99.9999999 

38.3612518  99.9999999  99-9999999 

35.3612518  279.999999  99.9999999 

WIND  LANES 

35.1779172  279.999999  72.7169939 

ship  lanes 

95.2255554  73.9999999  1.99999999 

16.5969849  279.999999  62.2442932 

LIST  or  HORIZONTAL  NOISE  ICUCLS  AND  T*€IR  ANGLES 
N  ,  DEGREES,  P9'-'U*U«3/'<KZ-DCQ> 

1  0.9999999996*99  35.9227242 

2  4.63364931  36.9226627 

3  45.9999999  32.3612612 

4  72.6999999  46.9326224 

5  73.5999999  32.3212619 

6  136.999999  36.2222627 

7  152.329692  36.9227222 

8  129.999999  36.9227222 

9  297.619492  36.9222227 

19  225.999999  32.3212619 

11  239.641922  49.9242627 

12  236.277263  49.2244674 

13  394.122279  49.9649267 

14  399.367219  32.3612619 

15  316.999999  36.9229127 

16  365.426399  36.9227249 

ONNI DIRECTIONAL  NOISE  D8//UPAM2/M2 

63.5114288  16.9999999  KNOTS 


49.9244674 


32.3612619 

36.9226627 

36.9227242 


01/11/88  NOISE  ROSE 

DUNES  (2.31  MINO  SPEED  =  15.0  KN  LEVELS  DB//UPRM2/ IHI-DECJ 


FREQUENCY 

504.0 

OMNI  NOISE 

DB//UPA«»2/HE 

63.5 


l  HZ ) 


ftKimmsis: 


t  to  database  files 


0  FILE  IN  OATA  BASES?  (V/N> 

Th€  UlHt  SPEEDS  K8UESTED  WILL  BE  OUTPUT 

CAL,  UNIT  9  •  OMNI ,  UNIT  19  •  HOftIZ.  UNIT  It  •  39 
EC  U0N3  (FO*  ALL).  UO  (FOB  ONLV  2)  V  OHJ 
CTIONAL  NOISE  DB"UPA4*2/MZ-S* 

1  S7. 5717921 

2  57.5717921 

3  93.5393958 

4  63.5393958 

5  74.9989981 

9  74.9889981 

7  79.5652826 

8  79.5552826 

9  74.9889981 

19  74.9889981 

It  63.5427246 

12  63.5427246 

13  58.4797974 

14  58.4797974 

29.9999999  77.2177429 

IRCCTIOHAL  NOISE  DB"UPA*42/MZ-S* 


1 

2 

3 

4 

s 

6 

7 

8 
9 

19 

11 

12 

13 

14 

26.9999999 


57.4881854 

57.4681854 

63.9383885 

63.9388885 

74.9993169 

74.9993169 

71.6599731 

71.6599731 

74.9993169 

74.9998169 

63.94C68S8 

63.9468858 

51.3763885 

58.3763885 
77.8751373 


IRECTIONAL  NOISE  Df//UPA8*2^MZ-5R 

1  57. 423 >3194 

2  57.4234194 

3  64.5195896 

4  64.5195896 

5  75.1233368 

6  75.1233368 

7  72.7567749 

8  72.7567749 

9  75.1233368 

19  76.1233368 

11  64.5297367 

12  64.5297367 

13  51.3327414 

14  SS.  3327484 

32.9999999  78.9315155 

IRECTIONAL  NOISE  D»/'''UPA*«'HZ-S* 

1  56.4258279 

2  52.4253279 

3  65.3389742 

4  25.3389749 

5  75.4237178 


Vertical  noise  levels 
at  14  angles  provided 
for  checking  3D  option 
data 


List  continues  for  all  frequencies  chosen  by 
DUNES  -  see  vertical  noise  values 


37 


WSRL-TN-34/89 


I. 11  Alternative  output  options  -  all  wind  speeds 


DUNES  2.3  OMNIDIRECTIONRL  NOISE 

PACIFIC  OCEAN  -  FEBRUARY  SHIP  ^ANE  PARAMETERS  MIND  uRNE  PARAMETERS 

LATITUDE  -30.0  DEGS 

M INOSPEEOS  5-40  KNOTS  TTPE  tEN  RANG  SHP3  LEVI  PEAK  BEAR  TTPE  LEN  NIOTh  RANK  BEAR  hSPO 

BOT .LOSS  GRADIENT*  1.0  SHIP*  0  100  I  165  60  73  SLPE  1000  25  100  270  15 

VERTICAL  BOT.LOSSi  0.0  LANE  1000  2CC  20  16S  60  270 

HlNOLANE  FRACTIONS! 

0.5  IN)  1.0  IS) 

HD  IHH200.DEPTMm.50  KH 

C3.  CR.  CB!  1527.  mBB.  1539.  QUAD.  RANGE  N.  5.  E.  N:  100.  100.  5000.  5000.  KH 


01/11/88  NOISE  ROSE 


DUNES  (2.31  HIND  5 PEED  =  *0.0  KN  LEVELS  D8//UPA«2/ IHZ-DEG) 


FREQUENCY  (HZ) 
504.0 

OMNI  NOISE 

D6//UPAnn2/HZ 

71.6 

qq7o!m  35.  KN 
3U69.  t  30.  KN 

67.6  25.  KN 

65.8  20.  KN 

63.5  15.  KN 

60.6  10.  KN 

58. 8  5.  KN 


180 
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COMMENTS : 

DATE 

FILE 

GDBMS  ADDRESS 


APPENDIX  II 

DUNES  2.3  RUNS  -  INPUT  WORKSHEET 

TE  ST 


OCEAN:  Pacific  [  V  ]  Indian  [ _ ]  Atlantic  [ 


RECEIVER  LATITUDE: 

-30 

LONGITUDE : 

(not 

used) 

RECEIVER  DEPTH  (M) : 

200 

„/ 

SEASON : 

February  [ 

V  1 

August  f  1 

SOUND  SPEEDS  (M/S): 

Default  [ 

]  or 

cs  1527  cr  1488  cb 

!S3^ 

BOTTOM  DEPTH  (KM) : 

45 

CLEAR  RANGE  OF  TRANSMISSION  BLOCKAGE  IN  4  QUADRANTS  (+-  45  DEG ) : 

N  0-1 

S  0-1 

E  S  W  S 

LOW  ANGLE  BOTTOM  LOSS 

SLOPE  FACTOR 

1 

(i 

.0  =  average) 

HIGH  ANGLE  BOTTOM  LOSS 

DB  SHIFT 

o 

(0 

dB  =  average) 

ANGLE  LIMIT  LOC  NOISE 

(Deg) 

or  CSUBB 

(M/S)  or  Default 

[  v  ] 

HIGH  LATITUDE  WIND  NOISE  FRACTIONS 

N  OS 

S  1 

AND  THEIR  BEARINGS  (DEG) 

N  O 

— 

S  180  (usually  0 

,180) 

WIND  LANES  1 

ENHANCED 

LENGTH 

WIDTH 

RANGE  BEARING 

WINDSPEED 

:  (Y=SLOPE ,N=STORM) 

!  (KM)  j 

(KM) 

1 

(KM)  |  (DEG)  RE  R| 

(KN) 

1  / 

iooo 

25 

* 

I0O  ;  270  i 

1  .  _ ' 

15  ’ 

i  i  :  i 

?  :  i  s  i  i  : 

j 


«  SHIP  CONTRIBUTIONS 


!  if  i 

ENHANCED  | 

i 

LENGTH 

RANGE  | 

BEARING  I 

#  SHIPS  | 

PEAK-LVL  | 

1 

FREQ  i 

!  i 

(Y/N) 

(KM)  i 

I  (KM)  | 

|  | 

(DEG)  RE  R 

1 

1 

(DB)  | 

-  -  \ 

(HZ)  | 

'  1' 

/ 

0 

'  100  1 

73  * 

1  ' 

165  ' 

6o  ' 

2 

N 

1000 

200  ; 

270 

20  : 

i65  ; 

- 1 

60 

i  i  ;  i  i  ! 

:  3i  ;  ;  !  i  i  i  ; 

fj  -  j  i  1  1  l  j 

i  i  i  i  1  i  t  i  i 

1  4|  III  1  1  I  1 

i  ;  i  I  •  i 

1  5!  II!  1  !  1  ! 

SELECTED  WIND  SPEED  1 5 

DATA 

OUTPUT  FOR 

DATABASES 

y  (y/n) 

WHICH  OF 

V  0  H  3  ? 

VOrt  3 
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DUNES  2.3  RUNS  -  INPUT  WORKSHEET 


COMMENTS : 

DATE 

FILE 

GDBMS  ADDRESS 


OCEAN : 


Pacific  [ _ ]  Indian  [ _ ]  Atlantic  [ _ ] 


RECEIVER  LATITUDE:  _  LONGITUDE:  _  (not  used) 

RECEIVER  DEPTH  (M) :  _ 

SEASON:  February  [ _ ]  August  [ _ ] 

SOUND  SPEEDS  (M/S):  Default  [  ]  or  CS  CR  CB 

BOTTOM  DEPTH  (KM):  _ 

CLEAR  RANGE  OF  TRANSMISSION  BLOCKAGE  IN  4  QUADRANTS  (+-  43  DEG): 

N  _  S  _  E  _  V  _ 

LOW  ANGLE  BOTTOM  LOSS  SLOPE  FACTOR  _  (1.0  =  average; 

HIGH  ANGLE  BOTTOM  LOSS  DB  SHIFT  _  (0  dB  =  average) 

ANGLE  LIMIT  LOC  NOISE  _  (Deg)  or  CSUBB  _ _  (M/S)  or  Default  [ _ 

HIGH  LATITUDE  WIND  NOISE  FRACTIONS  N  _ _  S  _ 

AND  THEIR  BEARINGS  (DEG)  N  _  S  _  (usually  0,180) 


#  WIND  LANES 


ENHANCED 

I  j  (Y=SL0PE .N-STORM)  ! 

LENGTH 

(KM) 

WIDTH  ' 
(KM)  | 

RANGE 

(KM) 

BEARING 
(DEG)  RE  R 

WINDS PEED 
(KN) 

'  1 


2 


!  3 

i  4i 

i - - - 

i  5| 

it  SHIP  CONTRIBUTIONS 


i  ENHANCED 

LENGTH 

RANGE  i 

BEARING 

,7  SHIPS  PEAK-LVL 

FREQ  , 

i  !  (Y/N)  j 

i _ l _ i 

(KM)  j 

(KM)  : 

(DEG)  RE  R 

i  (DB) 

(HZ) 

!  '  ' 


o 

3 


I  41 


SELECTED  WIND  SPEED 


DATA  OUTPUT  FOR  DATABASES 
WHICH  OF  V  0  H  3  ? 


(Y/N) 
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APPENDIX  III 

BOTTOM  LOSS  DATA  FOR  DUNES  2.3 

Finding  the  data  values  to  use  from  the  measured  or  theoretical  bottom  loss 
versus  grazing  angle  curve. 


Bottom 

Loss 

(dB) 


0  0  90 

c 

0  (Grazing  angle  (degrees)) 

Below  0^_  the  curve  is  approximated  by  a^  0  where  is  decibels/radian.  Above 

0  the  curve  is  a  constant  8- 
c 

III.  I  Input  data  requires  S,  where  a^  =  S  *  3  f 1 . 5  dB/rad. 

Default  is  S  =  1.0  For  a  high  slope  use  S  >  1. 

For  a  low  slope  use  S  <  1. 

111. 2  Input  data  requires  &2.  where  &  =  Bx  +  B2  (dB) 

and  8i  =  2.9  +  2  sin (3. 5  log  f  -  6.82) 

=  0.9  <  32  Hz 
=  3.88  >  500  Hz 

Default  is  B2  =  0  (dB)  For  high  loss  use  &2  >  0 

For  low  loss  use  8  2  <  0 
with  8  >  =  0.5  dB  always 

111. 3  Input  data  requires  0^,  the  breakpoint  from  slope  to  constant  value. 
It  is  found  by  either 

(i)  0^  specified,  suggested  values  are  the  angle  with  bottom  loss 

-3  dB  down  from  the  90°  value  or  angle  of 

intromission  when  C  /C,  ,  >  1 
o  bb 

(ii)  cos  1  (C  /C,,  ),  C,,  =  sound  speed  in  bottom  sediment 

o  bb  bb 

(iii)  default  cos  1  (C  / 1 .  l'-'-C.  ) ,  C,  =  sound  speed  at  bottom  C  is  water 

ODD  O 

sound  speed  at  receiver. 
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APPENDIX  IV 
NOTES  ON  IBM  PC  USEAGE 

Graphics  requires  the  equivalent  of  the  Color  Graphics  Adapter  and  a  graphics 
monitor.  More  capable  cards  and  monitors  will  have  the  minimum  requirement  of 
two  colours  on  a  640  by  200  resolution  screen. 

Graphics  dump  of  the  screen  goes  to  the  standard  IBM  PC  printer,  Epson  or 
similar  type.  Be  sure  to  use  the  GRAPHICS  command  before  running  DUNES  for 
the  first  time.  Add  it  to  your  ALTOEXE.BAT. 

Each  subroutine  is  compiled  using  MS-FORTRAN  3  and  above  versions.  Some 
subroutines  are  lumped  together  in  the  same  *.F0R  so  there  is  no  need  to  refer 
to  everyone  of  them. 

LINK  to  GRAFX+8087+FORTRAN  using  the  (5  option  for  the  OBJs  and  use  either 
DLINK.PLT  or  DLINK.NPL  for  plots  or  no  plots,  respectively. 

To  run  type  in 


DUNES , CON , TEST 1 . DAT 

to  use  the  data  file  TEST1.DAT  and  replies  go  to  the  terminal  CON. 

DUNES , CON , CON 

will  run  with  data  prompts  to  be  typed  in  as  you  go. 

The  DUNES.EXE  has  been  compressed  by  using  the  MS  supplied  program 
EXEPACK.EXE.  This  reduced  the  storage  from  172K  to  122K. 

The  8087  coprocessor  is  required. 

Output  to  database  files  will  generate  requests  for  file  names  for  UNIT 
numbers  8  to  11.  Give  answer  as  CON  or  TEST1 . OMN ,TEST1 . VER  accordingly. 
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APPENDIX  V 
SOURCE  LEVEL  UNITS 

For  source  levels  referenced  to  steradians  it  is  easy  to  show  that  these 
values  approximately  correspond  to  the  effective  source  level  that  would  be 
received  1  m  from  an  area  1  m  by  1  m.  First  we  examine  two  cases. 

(1)  The  vertical  angle  subtended  by  a  circular  plane  area  of  1  m2  at  a 
point  1  m  away  is  58.86°  (or  29.43°  half  angle).  This  corresponds  to  a 
solid  angle  of  0.811  sr.  The  ratio  of  the  plane  area  to  the  area  on  a 
sphere  centred  at  the  apex  of  the  angle  is  1:1.0690  or  -0.29  dB  re  1  sr. 

(2)  The  area  on  a  sphere  of  radius  1  m  centred  at  the  apex  by  the  angle  of 
1  sr  is  1  m* .  The  vertical  angle  subtended  is  65.54  (or  32.77°  half 
angle).  This  arrangement  corresponds  to  a  circular  area  which  is  0.84  m 
away.  The  ratio  of  the  plane  area  to  the  area  on  a  sphere  subtended  by 
1  sr  is  0.9204:1  or  -0.36  dB  re  1  sr. 

These  two  cases  bound  an  estimated  -0.3  dB  error  when  equating  the  two 
definitions  of  source  level.  This  value  could  be  subtracted  from 
equation  (1). 

The  vertical  linear  arrays  used  to  measure  the  vertical  noise  levels  had  a 
higher  end  fire  resolution  than  60°  due  to  more  than  7  hydrophones  or  adaptive 
beamforming.  Therefore  the  influence  of  the  source  directionality  over  the 
range  of  angles  between  ±  60°  and  ±  90°  can  be  shown  to  be  even  smaller  in 
determining  the  source  level. 
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Figure  1.  Transmission  types  to  receiver  at  R 


Sound  speeds: 

Ce  at  bottom 
Co  at  receiver 
Cs  at  surface 
CBB  in  subbottom 


Figure  2.  Vertical  angle  assignments  for  transmission  types 
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Frequency  (Hz) 

-  DUNES  algorithm  - Burgess  &  Kewley  measurement 


Figure  3(a).  Surface  source  level  of  wind  generated  noise  _  DUNES 

algorithm,  -  -  -  Burgess  and  Kewley  measurement 


to  20  30 

Windspeed  (knots) 


Figure  3(b),  Variation  of  source  depth  (m)  with  wind  speed  an 
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Frequency  (Hz) 


150 


Figure  5.  Dipole  source  decoupling  factor  versus  frequency 
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DIRECT  PATH  COMPONENT 
(30°  -  90°) 

Adjust  for  bottom  loss  by 
adding  top  panel  correction 


BOTTOM  BOUNCE  COMPONENT 
(12°-  30°) 

Adjust  for  waterdepth  H(m) 
and  bottom  loss  scalar  (L) 
by  adding: 

10 109  &L) 

RSR  COMPONENT 
(10°-  12°) 

Adjust  for  depth  excess  D(m) 
by  adding: 

10  log  (D/1000) 


RR  COMPONENT 
(±10°) 


Adjust  for  slope  windspeed 
W(kt)  and  range  R(km) 
by  adding: 


10  log 


100 


Figure  11.  Generalised  wind  noise  spectra 
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Figure  13(a).  Source  power  spectrum  level  of  a  600  ft  bulk  carrier  measured 
at  a  depression  angle  of  20° 


Figure  13(b).  Envelope  of  continuous  radiation  from  a  bulk  carrier  observed 
at  a  29°  depression  angle 
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